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Abstract 

The  High  Energy  Laser  End-to-End  Operational  Simulation  (HELEEOS)  off-axis 
scattering  algorithm  is  designed  to  predict  the  irradiance  that  will  be  detected  at 
a  given  off-axis  location  due  to  atmospheric  scattering  of  a  high-energy  laser.  The 
HELEEOS  system  models  the  propagation  of  the  laser  through  the  atmosphere,  ac¬ 
counting  for  such  effects  as  turbulence,  thermal  blooming,  and  atmospheric  absorp¬ 
tion.  The  HELEEOS  off-axis  scattering  algorithm  uses  the  scattering  phase  functions 
of  the  Mie  scattering  models  to  predict  the  amount  of  radiation  that  will  be  scattered 
toward  a  particular  observation  location  from  each  point  along  the  beam  path,  and 
the  total  irradiance  that  will  be  received  at  that  location.  Algorithm  outputs  were 
compared  with  data  from  a  laser  test  conducted  at  Wright-Patterson  Air  Force  Base 
using  1.53-1.59  /im  communications  lasers  operating  at  0.6  W  of  power.  The  off-axis 
irradiance  was  measured  using  an  Indium  Gallium  Arsenide  (InGaAs)  camera.  A 
software  application  was  developed  to  assist  in  camera  calibration  and  the  analysis  of 
the  collected  images.  The  application  was  used  to  reduce  the  data  from  each  image 
to  a  single  irradiance  value  which  could  be  compared  with  HELEEOS  predictions. 
Preliminary  results  show  an  agreement  within  1-2  orders  of  magnitude  between  the 
HELEEOS  algorithm  and  the  measured  off-axis  irradiance,  although  the  potential 
exists  for  improving  this  result  through  more  detailed  analysis  of  the  same  data  set. 
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EXPERIMENTAL  VALIDATION  TECHNIQUES  FOR  THE  HELEEOS 
OFF-AXIS  LASER  PROPAGATION  MODEL 


I.  Introduction 

The  purpose  of  this  thesis  is  to  develop  techniques  for  experimental  validating  the 
off-axis  propagation  algorithm  within  the  High  Energy  Laser  End-to-End  Operational 
Simulation  (HELEEOS)  model.  HELEEOS  is  a  software  package  developed  at  the  Air 
Force  Institute  of  Technology  to  simulate  laser  propagation  through  the  atmosphere. 
Its  core  purpose  is  to  predict  the  properties  of  a  high  energy  laser  beam  at  every  point 
along  its  path,  taking  into  account  its  interaction  with  the  atmosphere  [3].  In  2006, 
Scott  Belton  extended  the  model  to  simulate  off-axis  scattering,  giving  HELEEOS  the 
ability  to  determine  the  irradiance  received  from  laser  at  some  off-axis  location  [3]. 
Although  HELEEOS’s  basic  functionality  has  been  experimentally  validated,  the  off- 
axis  scattering  capability  has  not  previously  been  validated.  Tests  with  an  infrared 
camera  and  1-watt  infrared  lasers  at  Wright-Patterson  Air  Force  Base  provide  a  means 
to  determine  the  accuracy  of  the  off-axis  scattering  model.  HELEEOS  was  designed 
with  high-energy  laser  weapons  in  mind,  but  these  tests  provide  useful  data  for  veri¬ 
fying  the  off-axis  scattering  algorithm,  as  well  as  to  determine  the  viability  of  using 
the  same  model  for  lower  power  communication  lasers.  Additional  tests,  conducted  in 
December  2009  at  Dahlgren,  VA,  using  a  more  powerful  20-kilowatt  laser,  will  provide 
a  means  to  test  the  model  against  a  high  energy  laser. 

It  is  assumed  that  the  reader  has  a  basic  understanding  of  laser  propagation  and 
atmospheric  properties. 


1 


1.1  B  ackgr  ound 


Lasers  are  operated  in  the  atmosphere  for  a  variety  of  military  applications,  in¬ 
cluding  illumination  of  targets  for  bombs  and  missiles,  communication,  and  most 
recently  as  experimental  directed-energy  weapons.  Currently  the  Air  Force  is  build¬ 
ing  two  experimental  laser  weapons,  the  Airborne  Laser  (ABL)  and  the  Advanced 
Tactical  Laser  (ATL).  Both  use  a  Chemical  Oxygen  Iodine  Laser  (COIL).  The  ATL 
is  carried  inside  a  modified  C-130  aircraft  and  designed  to  damage  ground  targets 
along  slant  paths  of  up  to  10  km.  The  ABL,  on  the  other  hand,  is  intended  for  use 
against  theater  ballistic  missiles  such  as  the  SCUD,  and  would  be  operated  from  a 
modified  Boeing  747  aircraft.  The  ABL  operates  at  power  levels  on  the  order  of  one 
megawatt,  while  the  ATL  has  power  outputs  on  the  order  of  kilowatts  [2,  1], 

As  high-powered  lasers  become  more  common,  U.S.  government  agencies  will  re¬ 
quire  information  about  lasers  used  by  our  adversaries.  Atmospheric  scattering  makes 
it  possible  to  observe  a  laser  beam  from  an  off-axis  point  (i.e.  not  in  the  immediate 
path  of  the  beam),  opening  the  possibility  of  characterizing  the  laser  by  remote  ob¬ 
servation.  Molecules  and  aerosol  particles  cause  light  to  be  scattered  out  of  the  beam 
in  all  directions.  If  the  location  of  the  laser,  its  target,  and  the  observer’s  location  are 
known,  with  the  atmospheric  environment  sufficiently  well  characterized,  the  power  of 
the  laser  can  be  determined  from  the  irradiance  of  scattered  light  seen  by  an  off-axis 
observer.  It  may  also  be  possible  to  use  the  same  technique  on  some  communications 
lasers.  A  computer  program,  written  by  Belton  and  since  integrated  into  HELEEOS, 
provides  the  computational  tools  to  determine  off-axis  scattered  irradiance  for  a  laser 
of  known  properties  [3] . 

HELEEOS  predicts  laser  propagation  by  an  empirical  method,  taking  into  account 
a  variety  of  factors  including  the  atmospheric  environment,  and  the  laser’s  ability  to 
compensate  for  atmospheric  turbulence.  The  model  supports  scattering  calculations 
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for  electromagnetic  beams  operating  between  355  nm  and  8.6  m.  HELEEOS  provides 
the  intensity  of  light  in  all  directions  at  each  point  along  the  beam;  Belton’s  model 
uses  this  to  determine  the  irradiance  received  at  a  particular  off-axis  point,  taking 
into  account  extinction  of  light  traveling  from  the  beam  to  the  observer  [3] . 

While  Belton’s  model  predicts  the  off-axis  scattered  irradiance  for  a  laser  whose 
properties  are  known,  a  new  feature  of  HELEEOS  due  to  be  added  by  Eric  Magee  and 
William  Gruner  is  the  ability  to  determine  the  power  of  a  laser  based  on  measurements 
of  off-axis  scattering.  This  feature  will  be  implemented  using  a  Gaussian  process 
regression  model  to  map  off-axis  irradiances  to  laser  power  outputs.  The  Gaussian 
process  technique  also  provides  the  ability  to  determine  a  probability  distribution  to 
indicate  the  degree  of  uncertainty  in  the  prediction  [11]. 

1.2  Motivation 

The  capability  to  determine  laser  power  from  backscatter  would  provide  a  means 
to  determine  basic  characteristics  of  noncooperative  laser  systems.  In  the  case  of 
communications  lasers,  it  may  also  be  possible  to  collect  some  portion  of  the  data 
being  transmitted.  Thus  the  technique  would  provide  an  important  intelligence  gath¬ 
ering  capability,  particularly  if  laser  weapons  entered  widespread  use  among  U.S. 
adversaries. 

A  software  model  already  exists  as  a  component  of  the  HELEEOS  computer  pro¬ 
gram.  The  core  functionality  of  the  HELEEOS  program,  that  is,  to  characterize  laser 
beam  propagation  for  a  known  laser  system,  was  validated  in  previous  studies.  This 
current  study  provides  an  initial  validation  for  the  off-axis  scattering  model  developed 
by  Belton,  and  develops  techniques  which  can  be  used  for  future  testing. 
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1.3  Problem  Statement 


The  primary  goal  of  this  research  is  to  develop  techniques  for  experimental  val¬ 
idation  of  the  HELEEOS  off-axis  scattering  model,  and  to  provide  a  preliminary 
assessment  of  its  accuracy.  A  secondary  goal  is  to  determine  whether  the  techniques 
used  to  measure  off-axis  irradiance  are  effective  for  communications  lasers. 

1.4  Methodology 

Tests  conducted  in  July  2009  by  the  Air  Force  Research  Laboratory  (AFRL)  at 
Wright-Patterson  Air  Force  Base  provided  the  first  set  of  validation  data.  The  lasers 
used  in  the  test  operated  at  1.533,  1.589,  and  1.595  pm,  within  the  range  supported 
by  HELEEOS,  and  with  power  levels  varying  from  0.6  to  5.7  W.  This  comparatively 
low  power  level  may  produce  results  that  differ  somewhat  from  megawatt-class  COIL 
lasers  of  interest,  but  testing  against  low  powers  opens  up  the  possibility  of  charac¬ 
terizing  noncooperative  communications  lasers  as  well  as  laser  weapons.  The  more 
powerful  weapons  lasers  would  likely  be  easier  to  detect  and  measure  than  these 
communications  lasers. 

Additional  tests  were  conducted  in  December  2009  using  a  1.07  pm  laser  in 
Dahlgren,  VA  [11],  These  tests  used  a  20- kilowatt  laser  intended  for  use  as  a  weapon, 
so  they  help  fulfill  the  need  to  demonstrate  the  off-axis  scattering  technique  against 
a  high-energy  laser. 

Irradiance  measurements  at  the  July  2009  test  were  made  with  an  infrared  camera 
using  an  Indium-Gallium- Arsenide  (InGaAs)  array.  Although  the  HELEEOS  off-axis 
scattering  module  was  designed  with  a  single-element  detector  in  mind,  the  camera 
provided  a  more  flexible  method  of  data  collection.  By  analyzing  the  images  collected 
it  is  possible  to  clearly  distinguish  background  pixels  from  laser  beam  pixels,  and 
the  background  irradiance  can  easily  be  subtracted  out  to  find  the  irradiance  of  the 
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scattered  laser  light  [11],  Furthermore,  decisions  could  be  made  after  the  data  was 
collected  about  which  part  of  the  camera’s  held  of  view  should  be  analyzed,  which 
might  not  have  been  possible  using  a  single-element  detector. 

Measurements  collected  during  the  December  2009  test  used  a  CCD  camera.  Al¬ 
though  CCD’s  are  designed  to  detect  visible  light,  they  have  adequate  sensitivity  at 
the  1.07  jirn  laser  wavelength  to  detect  the  scattered  laser  light. 

1.5  Overview 

HELEEOS  simulates  a  number  of  mechanisms  by  which  the  atmosphere  affects 
laser  beam  propagation.  All  of  these  are  described  in  Chapter  II.  The  most  important 
of  these  is  atmospheric  scattering,  which  will  be  discussed  in  greater  detail.  The 
scattered  light  from  the  optical  scattering  was  collected  by  an  InGaAs  camera,  which 
was  calibrated  against  a  blackbody.  Chapter  III  discusses  the  radiometry  theory 
required  to  calibrate  the  camera  and  apply  the  calibration  to  field  images,  while  the 
results  of  the  calibration  are  described  in  Chapter  IV.  The  calibration  was  applied  to 
56  test  cases,  which  are  compared  against  HELEEOS  runs  in  Chapter  V.  For  those 
interested,  a  detailed  summary  of  the  input  data  and  results  for  all  the  test  cases  can 
be  found  in  Appendices  B  and  C. 
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II.  Literature  Review 


This  chapter  summarizes  some  of  the  fundamental  concepts  relating  to  atmo¬ 
spheric  laser  propagation  and  scattering  detection  using  infrared  cameras.  The  Gaus¬ 
sian  process  regression  to  be  used  for  correlating  scattered  irradiance  with  laser  power 
is  also  described.  Finally,  a  discussion  of  communications  lasers  is  included  as  a  pre¬ 
lude  for  later  discussion  on  the  viability  of  using  HELEEOS  against  such  lasers. 

2.1  Atmospheric  Interaction 

Lasers  interact  with  the  atmosphere  in  a  number  of  ways.  Molecules  in  the  air 
(primarily  C02,  N2,  H20,  and  02 )  may  both  scatter  and  absorb  radiation.  Similarly, 
aerosols  may  contribute  both  scattering  and  absorption.  Finally,  optical  turbulence 
(variances  in  the  index  of  refraction  of  the  air)  can  perturb  the  beam  position,  expand¬ 
ing  it  over  a  time  average.  Turbulence  can  also  cause  localized  irradiance  fluctuations 
in  the  beam  [22,  18]. 

For  this  thesis,  the  interactions  of  primary  concern  are  absorption  and  scattering. 
Scattering  along  the  beam  is  necessary  for  light  to  reach  any  off-axis  detector,  but 
any  additional  scattering  along  the  path  between  the  beam  and  the  detector  will 
contribute  to  extinction,  decreasing  the  observed  irradiance.  Absorption  decreases 
the  irradiance  both  of  the  beam  itself  and  of  the  off-axis  scattering. 

Regardless  of  the  mechanism,  a  laser  beam’s  power  decays  exponentially  with 
distance  as  long  as  the  atmospheric  properties  remain  constant.  The  decay  rate  7  is 
given  by 


7  —  am  +  Pm  +  <Ti  +  Pai  (1) 

where  arn  and  prn  denote  absorption  and  scattering,  respectively,  by  molecules,  and 
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aa  and  /3a  denote  absorption  and  scattering  by  aerosol  particles  [18].  The  degree  of 
attenuation  at  a  particular  distance  is  then  given  by  Beer’s  Law: 


r  = 


/n 


=  e“7L 


(2) 


where  r  is  the  transmittance,  7  is  the  attenuation  coefficient,  and  2  is  the  distance 
from  the  source  [18].  If  the  atmospheric  properties  vary  along  the  beam  path,  the 
exponent  —7 z  is  replaced  by  the  integral 


-  [  7 (z)dz.  (3) 

Jo 

The  absorption  and  scattering  coefficients  depend  on  the  concentrations  of  the 
contributing  particles,  and  may  be  expressed  as  the  product  of  a  cross  section  with 
the  particle  concentration  [18].  The  cross  sections  are  given  as  cra  for  absorption  and 
as  for  scattering,  and  the  absorption  and  scattering  coefficients  are  given  by 

a  =  a aN a  (4) 

and 

(3  =  asNs,  (5) 


where  Na  and  Ns  are  the  concentrations  of  absorbers  and  scatterers,  respectively. 

2.1.1  Absorption. 

Each  molecular  species  in  the  atmosphere  has  a  characteristic  absorption  spectrum 
which  contributes  to  the  total  attenuation.  Of  the  molecular  species  present  in  the  at¬ 
mosphere,  H20  and  C02  contribute  the  most  to  absorption.  These  molecules  absorb 
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incident  light  by  changing  between  various  rotational  and  vibrational  modes  [18]. 

Many  high  energy  lasers  are  designed  to  operate  in  “clean- window”  bands  for 
which  molecular  absorption  is  minimized.  For  these  lasers,  scattering  is  the  dominant 
attenuation  mechanism  [10]. 

2.1.2  Scattering. 

Atmospheric  scattering  is  commonly  treated  using  two  models:  Rayleigh  scatter¬ 
ing  and  Mie  scattering.  Rayleigh  scattering  models  scattering  by  particles,  primar¬ 
ily  molecules,  which  are  smaller  than  the  size  of  the  wavelength  in  question.  Mie 
scattering  models  larger  particles  for  which  the  shape  of  the  scattering  particle  is  a 
factor  [18]. 

Rayleigh  scattering.  Rayleigh  scattering  occurs  when  the  electric  held  of 
the  laser  beam  displaces  bound  electrons  in  molecules  of  the  air.  This  causes  a 
dipole  moment  which  oscillates  at  the  same  frequency  as  the  incident  held.  For  light 
scattered  in  a  direction  D,  the  scattered  intensity  has  the  proportionality 

I  oc  u4(l  —  sin2(6l)  cos2  0),  (6) 

where  6  is  the  angle  between  Q  and  the  direction  of  propagation,  and  u  is  the  frequency 
of  the  light.  0  is  the  angle  between  the  electric  held  of  the  polarized  light  and  the 
component  of  Q  in  a  plane  perpendicular  to  the  propagation  direction  [14] .  Dropping 
the  frequency  dependence  gives  us  the  scattering  phase  function 

p(0, 0)  =  cn4(l  —  sin2(d)  cos2  0).  (7) 

This  scattering  prohle,  also  known  as  the  scattering  phase  function,  is  shown  in 
Figure  1. 
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Figure  1.  Scattering  phase  function  for  Rayleigh  scattering  in  directions  6  relative 
to  the  propagation  direction.  Results  are  shown  for  light  polarized  perpendicular, 
parallel,  and  unpolarized  relative  to  the  scattering  direction  (based  on  Figure  12.2  and 
Equation  12.9  of  [14]). 


Although  the  Rayleigh  scattering  model  accounts  for  polarization  effects,  HELEEOS 
threats  the  light  as  being  unpolarized  during  propagation,  except  for  considering  po¬ 
larization  that  is  induced  by  the  scattering  process  itself  [11].  Furthermore,  no  polar¬ 
ization  information  was  available  for  the  lasers  used  in  the  test  data  for  this  project. 
Studies  previously  conducted  by  the  AFIT  Center  for  Directed  Energy  show  that 
polarization  has  little  effect  on  HELEEOS  outputs  for  lasers  below  12  km  of  altitude, 
as  shown  in  Figure  2.  This  is  because  the  dominant  scattering  mechanism  at  lower 
altitudes  is  aerosol  scattering,  which  does  not  induce  polarization  [11], 

The  frequency  dependence  of  Rayleigh  scattering  means  that  high-frequency  light 
is  scattered  most  strongly.  This  effect  is  responsible  for  the  blue  color  of  the  sky,  since 
blue  light  is  scattered  more  strongly  than  red.  The  frequency  dependence  also  means 
that  Rayleigh  scattering  plays  a  minimal  role  for  wavelengths  greater  than  1  jim. 


9 


Figure  2.  Effect  of  polarization  on  HELEEOS  off-axis  irradiance  outputs.  Although 
polarization  can  make  large  differences  at  high  altitudes,  it  has  little  effect  below  1,200 
km. 


Mie  scattering.  Aerosol  scattering  is  the  most  important  factor  for  lasers 
operating  in  the  troposphere  (which  is  the  case  for  all  the  planned  validation  tests)  [10] . 
Aerosols  are  small  particles  in  the  air  which  can  include  dust,  salt  grains,  and  water 
droplets.  Rayleigh  scattering  ignores  the  geometry  of  the  scatterers  and  treats  them 
as  point  particles,  but  this  assumption  is  invalid  for  scatterers  whose  size  is  on  the 
order  of  the  laser  wavelength  or  larger.  Therefore  the  Mie  model  is  instead  used  for 
these  particles. 

The  attenuation  cross-section  a  for  Mie  scattering  is  given  by 

a  =  Kira 2,  (8) 

where  K  is  called  the  attenuation  factor  and  a  is  the  radius  of  the  scattering  particle. 
In  general  K  may  include  both  an  attenuation  term  and  a  scattering  term,  and 
depends  on  both  wavelength  and  the  radius  a: 
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K{a,  A)  Kscaf(^d^  A)  -\-  A). 


(9) 


Kscat  and  Kai)S  are  also  called  absorption  and  scattering  efficiencies. 

The  presence  of  both  absorption  and  scattering  terms  in  K  presents  a  challenge, 
particularly  in  the  infrared  where  the  molecular  absorption  coefficients  can  be  sig¬ 
nificant.  This  is  further  compficated  by  the  fact  that  the  particle  sizes  may  vary, 
and  normally  there  are  multiple  particle  sizes  present.  In  order  to  determine  the  Mie 
attenuation  coefficient  it  is  then  necessary  to  obtain  a  particle  size  distribution  and 
integrate  that  distribution  over  the  various  radii  [18]: 

/3(A)  =  7T  /  N  (a)  K scat  a2  da.  (10) 

J a\ 

Normally  the  Mie  computations  are  performed  using  a  nondimensional  size  pa¬ 
rameter 


27ra 

so  that  the  results  of  Mie  scattering  computations  can  be  compared  across  wavelengths 
and  the  results  applied  more  generally.  In  the  case  of  x  «  1,  the  Mie  model  predicts 
Rayleigh  scattering  behavior,  in  which  forward  and  back  scattering  are  roughly  equal. 
At  larger  size  parameters,  forward  scattering  is  stronger,  as  shown  in  Figure  3. 
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Figure  3.  Mie  scattering  phase  functions  for  various  size  parameters  (based  on  Fig¬ 
ure  12.7  of  [14]).  The  radial  axis  is  a  logarithmic  scale,  and  the  angles  are  relative  to 
the  direction  of  propagation.  Features  of  note  include  the  peaks  at  0°  and  180°,  corre¬ 
sponding  to  forward  and  back  scattering,  respectively,  as  well  as  the  fact  that  all  the 
phase  functions  reach  their  minimum  values  somewhere  around  90°.  Finally,  it  should 
be  noted  that  the  scattered  intensity  increases  with  the  size  of  the  scattering  particles, 
denoted  by  the  size  parameter  x. 


The  Mie  model  makes  the  simplifying  assumption  that  all  the  scatterers  are  spher¬ 
ical  in  shape.  The  spherical  particle  assumption  can  introduce  significant  errors  at 
short  wavelengths,  but  according  to  Yang  et  al.  the  effect  of  particle  nonspheric¬ 
ity  is  insignificant  at  thermal  infrared  wavelengths  (8-12  / im )  [21].  Although  the 
wavelengths  of  interest  in  this  study  are  shorter  than  that,  the  assumption  may  still 
be  relevant  because  random  variations  of  the  orientations  and  sizes  of  the  particles 
cause  the  errors  to  average  out  [11].  The  HELEEOS  model  therefore  uses  Mie  scat¬ 
tering  exclusively  to  simulate  atmospheric  scattering.  The  phase  function  used  for 
HELEEOS  runs  in  this  thesis  is  shown  in  Figure  4;  this  phase  function  is  computed 
within  HELEEOS  using  an  implementation  of  the  Wiscombe  Mie  algorithm  [10],  a 
technique  published  in  1980  which  provides  superior  performance  compared  to  previ¬ 
ous  methods  of  computing  Mie  scattering  phase  functions  [19]. 
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Figure  4.  Mie  scattering  phase  functions  used  in  HELEEOS.  Both  phase  functions  are 
computed  from  the  Wiscombe  Mie  module,  but  with  different  size  parameters.  Note 
that  aerosol  scattering  dominates;  this  is  because  the  Mie  model  predicts  stronger 
scattering  for  larger  particles. 


2.1.3  Optical  turbulence. 

Optical  turbulence  is  caused  by  random  localized  temperature  variations  in  the  air, 
which  affect  the  index  of  refraction  and  thus  distort  the  shape  of  the  beam,  causing 
it  to  be  wider  over  a  time  average  than  it  would  be  in  a  vacuum.  Small  pockets  of 
turbulence  (smaller  than  the  beam  diameter)  can  also  cause  intensity  fluctuations,  or 
scintillation,  in  the  beam  [18].  Optical  turbulence  is  an  important  factor  influencing 
laser  beam  quality  and  spread,  but  it  is  of  less  importance  for  off-axis  scattering.  The 
broadened  beam  will  impose  lower  irradiances  on  the  scattering  particles  in  the  air, 
but  no  loss  will  result  when  one  looks  at  the  entire  beam  as  a  whole. 
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2.2  Infrared  detection 


2.2.1  InGaAs  detector  properties. 

The  data  used  in  this  thesis  were  collected  using  an  infrared  camera  containing  an 
indium-gallium-arsenide  (InGaAs)  focal  plane  array.  An  InGaAs  array  benefits  from 
a  low  dark  current  and  high  quantum  efficiency.  In  addition  to  cameras,  InGaAs 
detectors  are  often  used  for  fiber-optic  communication. 

InGaAs  arrays  are  popular  infrared  detectors  because  of  their  wide  bandwith. 
InGaAs  detectors  have  a  response  range  from  0.5  to  1.7  /im,  with  the  strongest 
response  between  1.2  and  1.5  firn,  as  shown  in  Figure  5  [9].  This  overall  response 
range  overlaps  well  with  the  1.5-1. 6  /im  wavelength  range  of  the  lasers  in  question. 


Figure  5.  Responsivity  spectrum  for  a  typical  InGaAs  detector  at  25°C  with  a  5V  bias. 
Based  on  Figure  7.55  in  [9]. 


2.2.2  CCD  detector  properties. 

The  data  from  the  December  2009  test  were  collected  using  a  charge-coupled  device 
(CCD),  which  employs  an  array  of  silicon  photodiodes.  Silicon  detects  radiation 
between  0.2  and  1.15  /ira,  with  the  best  response  between  0.6  and  0.8  /jm,  as  shown 
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in  Figure  6.  Because  this  overlaps  well  with  the  the  response  range  of  human  eyesight, 
the  silicon  photodiode  is  the  most  commonly  used  type  of  photovoltaic  detector  [9]. 
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Figure  6.  Relative  response  spectra  for  CCD  detectors.  Based  on  Figure  7.30  in  [9]. 

Silicon’s  spectral  response  would  make  it  unsuitable  for  the  July  2009  test,  but 
silicon  does  respond  at  the  shorter  1.05  /im  wavelength  of  the  laser  used  in  the  De¬ 
cember  2009  test.  Therefore,  although  it  is  less  suitable  for  infrared  detection  than 
an  InGaAs  array,  the  CCD  is  capable  of  detecting  the  radiation  of  interest  for  in  the 
case  of  the  December  test. 

2.3  Gaussian  process  regression 

Gaussian  process  regression  is  a  technique  for  “machine  learning,”  a  process  of 
mapping  some  set  of  input  parameters  to  an  output  parameter.  A  machine  learning 
process  takes  a  set  of  input  data  with  corresponding  output  data  (the  “training” 
data  set)  and  produces  a  numerical  function  which  maps  the  inputs  to  the  outputs, 
attempting  to  smooth  out  noise  such  that  the  function  will  work  for  other  inputs 
besides  those  used  in  the  training  set  [15]. 

A  machine  learning  process  should  fit  the  training  data  with  reasonable  accuracy, 
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but  not  so  well  that  noise  in  the  training  data  influences  the  training  data  too  much. 
If  the  fit  is  too  good,  the  model  may  be  unstable,  meaning  that  it  works  excellently 
for  the  training  data  but  poorly  for  any  new  (or  “test”)  data  with  which  it  is  used. 

A  simple  example  of  machine  learning  is  linear  regression,  in  which  the  mapping 
function  is  simply  the  sum  of  one  or  more  known  functions,  which  are  scaled  by 
coefficients.  This  approach  is  relatively  simple  and  its  output  easy  to  interpret.  It 
affords  the  implementor  tight  control  over  basis  functions,  which  is  desirable  when 
the  physics  of  the  system  to  be  modeled  is  well  known.  The  linear  model  cannot, 
however,  model  systems  involving  complex  interactions  between  input  parameters. 
One  can  always  improve  the  quality  of  the  fit  to  the  input  data  by  increasing  the 
number  and  variety  of  basis  functions,  but  this  can  lead  to  over- fitting  if  one  is  not 
careful. 

More  complex  systems  that  cannot  be  modeled  by  linear  regression  are  sometimes 
modeled  using  neural  networks.  Neural  networks  are  more  flexible  than  the  linear 
models,  but  they  are  in  some  respects  too  flexible.  Neural  networks  offer  a  wide  range 
of  choices  for  the  implementor  in  terms  of  design  parameters,  and  no  straightforward 
framework  yet  exists  to  aid  in  such  decisions  [15].  As  with  linear  models,  the  high 
degree  of  flexibility  can  lead  to  overfitting  if  one  is  not  careful  [12] . 

Gaussian  process  regression  offers  an  alternative  which  is  more  flexible  than  the 
linear  models  but  provides  a  rigorous  approach  to  avoiding  an  over-fit  of  the  training 
data.  It  can  work  with  a  large  (even  infinite)  set  of  basis  functions,  while  maintaining 
the  stability  of  the  model  by  computing  a  characteristic  length  scale.  The  character¬ 
istic  length  scale  is  a  measure  of  how  far  one  can  move  through  the  input  data  set 
before  seeing  a  significant  change  in  the  output  [15]. 

In  the  HELEEOS  off-axis  scattering  model,  Gaussian  process  regression  will  be 
used  to  create  a  model  that  predicts  the  laser  power  for  a  given  off-axis  observation 
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scenario.  A  training  data  set  will  initially  be  produced  by  running  the  existing  off-axis 
scattering  model  with  various  input  parameters.  After  more  experimental  data  are 
available,  they  will  be  used  to  form  a  new  training  data  set. 

2.4  Laser  communications 

Although  they  are  not  a  predominant  means  of  communication  today,  lasers  are 
of  interest  for  their  high  bandwidth,  security,  and  freedom  from  radio  interference 
difficulties  [13,  4],  The  line-of-sight  nature  of  laser  beams  means  that  data  carried 
through  them  will  not  interfere  with  radio  communications,  and  as  a  result  laser 
communications  are  not  currently  regulated  by  the  FCC  [4],  Since  nearly  all  the 
transmitted  power  is  directed  at  the  intended  receiver,  it  is  possible  to  obtain  a  desired 
signal  strength  at  the  receiver’s  end  using  a  much  smaller,  lighter,  and  less  powerful 
transmitter  [13].  For  military  applications,  the  line-of-sight  restriction  means  that 
laser  communications  are  more  secure  than  radio  communications  because  they  are 
more  difficult  to  intercept  [13]. 

The  particular  lasers  used  for  the  test  cases  in  this  study  were  designed  for  use 
in  laser  communications,  and  operated  at  power  levels  of  0.6  to  5.7  watts.  Thus  the 
data  already  demonstrates  the  possibility  of  detecting  the  use  of  laser  communications 
systems.  Depending  on  the  manner  in  which  laser  communications  were  transmitted, 
it  might  also  be  possible  to  decipher  some  of  the  transmitted  data.  Barring  that, 
software  such  as  HELEEOS  could  be  used  to  glean  information  on  the  communications 
system  itself,  such  as  its  operating  power.  Doing  so  requires  that  the  camera  observing 
the  laser  operations  be  carefully  calibrated.  The  following  chapter  describes  the  theory 
required  to  calibrate  the  camera,  and  to  use  the  calibration  result  to  compute  a  single 
ir radiance  value  for  an  image. 
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III.  Theory 


In  order  to  convert  the  raw  camera  data  into  irradiances  which  could  be  compared 
with  HELEEOS,  the  camera  had  to  be  calibrated  against  a  known  radiation  source, 
in  this  case  a  blackbody  whose  spectral  radiance  is  described  by  Planck  theory  Ra¬ 
diometric  theory  allows  this  radiance  to  be  converted  into  the  power  received  at  the 
camera’s  aperture  as  long  as  the  aperture  size,  the  distance  and  direction  from  the 
blackbody  to  the  camera,  and  the  camera’s  field  of  view  are  all  known. 

Once  the  calibration  has  been  applied  to  images  taken  in  the  field,  geometric 
conversions  must  be  made  in  order  to  determine  the  distance  along  the  laser  beam 
path  that  corresponds  to  each  pixel.  In  order  to  accomplish  this,  the  angle  between 
the  laser  beam  and  a  line  from  laser  source  to  the  observer,  referred  to  as  the  off-axis 
angle,  must  be  known.  Once  this  is  done,  the  irradiances  computed  from  the  camera 
can  then  be  compared  with  irradiance  values  from  HELEEOS  in  a  meaningful  way. 


3.1  Camera  Calibration 


In  order  to  convert  the  digital  outputs  from  the  camera  into  irradiances,  the 
camera  must  be  calibrated  against  a  known  radiation  source,  and  a  blackbody  was 
chosen  for  this  purpose.  The  radiance  of  the  blackbody  at  a  particular  temperature 
and  wavelength  is  given  by  the  Planck  radiation  law  as  [6] 


Lx  = 


2  he2 


(12) 


A"  exP  (wr)  ~  1 

The  form  of  the  Planck  formula  given  here  is  the  spectral  (wavelength  dependent) 
form.  Strictly  speaking,  photodetectors  such  as  an  InGaAs  array  respond  to  photons 
rather  than  energy.  That  is,  a  single  photon  will  liberate  one  electron  if  it  is  sufficiently 
energetic,  and  the  energy  of  individual  photons  is  immaterial  as  long  as  it  is  above 
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that  threshold.  Therefore,  it  would  be  more  correct  to  use  the  photon  count  form  of 
the  Planck  law,  rather  than  the  radiant  flux  form  used  here.  However,  for  a  narrow 
wavelength  range  the  two  equations  produce  essentially  the  same  result.  The  July 
2009  measurements  and  the  corresponding  calibration  were  conducted  using  a  filter 
with  sufficiently  narrow  band  that  the  energy  form  of  the  Planck  equation  may  be 
safely  substituted  for  the  photon  count  form.  The  energy  form  was  chosen  over  the 
photon  count  form  because  it  produces  outputs  in  the  form  of  an  irradiance  rather 
than  a  photon  flux,  so  that  the  units  are  the  same  as  the  output  from  HELEEOS. 

The  radiance  L  is  related  to  the  exitance  M  by 

M  =  -^T=f  Lcos6sdttd.  (13) 

oAs  J<nd 

For  a  blackbody,  the  radiance  L  is  independent  of  the  direction  of  Qd.  A  radiation 
source  having  this  property  is  called  a  Lambertian  radiator.  For  a  planar  Lambertian 
source,  the  relation  simplifies  to  [9] 


M 


f 
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^ d 
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Ultimately  the  goal  is  to  find  the  power  (j>  received  at  the  aperture,  which  is  related 
to  the  radiance  L  by 


L 


cP(p 

dAproj  dfly/ 


(15) 


where  cL4pro j  is  an  area  element  of  the  emitting  surface  projected  into  a  plane  perpen- 
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dicular  to  the  line  of  transmission.  eL4pr0j  =  dAs  cos  6S  where  dAs  is  the  area  of  the 

emitting  surface  and  0S  is  the  angle  between  the  line  of  transmission  and  the  normal 

to  the  emitting  surface,  and  dAlj  is  the  solid  angle  subtended  by  the  camera  aperture 
as  viewed  from  the  surface  element  in  question,  and  whose  value  is  given  by 

«  =  £>  (16) 

where  a  is  the  surface  area  of  the  aperture  and  r  is  the  distance  from  the  emitting 
surface  element  to  the  aperture.  Solving  equation  15  for  d2<f>  and  integrating  yields 

cj)  =  J  j  LdAprojdCl.  (17) 

In  the  calibration  tests,  the  blackbody  has  an  angular  extent  of  ±3°  from  the 
vantage  point  of  the  camera,  so  by  a  small  angle  approximation  dAproj  ~  dA.  [16] 
Similarly  the  solid  angle  fl  can  be  treated  with  a  constant  radius  r  =  600cm.  Thus 
equation  17  simplifies  to 


4>  =  LAVL  (18) 

in  the  small-angle  approximation. 

To  further  verify  that  this  small-angle  approximation  is  valid,  the  full  integral 
of  equation  17  can  be  performed  numerically.  For  the  square  planar  blackbody  in 
question,  this  becomes 
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and 


vdc2  +  y2  +  r 
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so  that  the  integral  to  be  evaluated  is 


(21) 
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Using  the  values  of  a  =  ^-7r  cm,  r  =  600  cm,  and  w  =  15.24  cm,  the  integral  22 
evaluates  to  4>  =  0.206624 L,  while  the  small-angle  approximation  of  equation  18 
evaluates  to  0  =  0.206757L.  Thus  the  small-angle  approximation  introduces  an  error 
of  about  0.06%.  Note  that  this  numeric  integral  still  uses  a  small-angle  approximation 
for  the  evaluation  of  U  inside  the  integral;  this  introduces  an  additional  error  but  this 
error  will  be  smaller  than  0.06%  since  the  aperture  of  the  telescope  is  smaller  than 
the  width  of  the  blackbody.  Therefore  the  total  error  caused  by  the  small  angle 
approximation  is  likely  to  be  at  most  0.1%. 

The  radiant  flux  (or  power)  is  multiplied  by  the  integration  time  (a  parameter 
controlled  by  the  camera  operator)  to  obtain  a  radiant  energy  Q.  This  can  then  be 
divided  by  the  total  number  of  pixels  in  the  array  to  obtain  a  radiant  energy  received 
per  pixel. 


3.2  Experiment  geometry  calculations 

In  analyzing  the  data,  it  will  be  necessary  to  relate  pixels  in  the  scene  to  physical 
locations  in  the  experiment.  In  order  to  do  this,  the  geometry  relevant  to  each  pixel 
is  considered  relative  to  the  source  and  observer  as  shown  in  Figure  7.  The  location 
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of  a  pixel  along  the  laser  beam  path  is  described  relative  to  the  source  in  the  form 
of  the  distance  BD,  and  relative  to  the  observer  in  terms  of  the  angle  o,  called  the 
observer  angle.  In  order  to  determine  the  distance  BD  from  the  angle  a,  or  vice  versa, 
the  off-axis  angle  b  and  the  distance  AB  must  be  known. 


Observer 


Figure  7.  Experiment  setup  with  a  triangle  shown  between  the  observer,  source,  and 
a  point  of  interest  along  the  laser  beam  path. 


3.2.1  Determining  the  off-axis  angle. 

Due  to  the  extreme  distances  involved,  it  is  seldom  possible  to  measure  the  off-axis 
angle  b  directly.  Instead,  the  angle  must  be  determined  from  the  distances  between 
source,  observer,  and  target.  The  angle  is  then  given  by  the  Law  of  Cosines  as 

(BC)2  +  (AB)2  -  (BC)2 
b~  2(AB)(BC) 

where  BC  is  the  distance  from  the  source  to  the  target,  AB  is  the  distance  from 
the  source  to  the  observer,  and  AC  is  the  distance  from  the  observer  to  the  target, 
relating  to  the  points  A,  B,  C,  and  D  in  Figure  7. 


3.2.2  Entering  observer  geometry  into  HELEEOS. 

The  distance  AB  and  the  angle  b  are  entered  into  the  observer  geometry  tab  of 
HELEEOS,  shown  in  Figure  8.  The  distance  AB  is  entered  into  the  “Initial  Distance 
from  Platform  (m)”  field.  In  the  scenarios  used  here,  the  azimuth  from  the  Platform 
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to  the  Target  was  set  to  90°,  so  the  off-axis  angle  b  is  not  entered  directly,  but  rather 
90  —  b  is  entered  into  the  Relative  Azimuth  from  Platform  field.  In  practice  this  data 
entry  process  was  automated  by  a  script. 


Figure  8.  HELEEOS  observer  geometry  settings.  The  settings  in  the  Initializations 
section  are  changed  according  to  the  location  of  the  observer. 

3.2.3  Locating  a  pixel  position  along  the  beam. 

The  HELEEOS  off-axis  scattering  model  requires  a  specific  point  or  range  of  points 
along  the  length  of  the  beam  at  which  to  compute  the  scattering  toward  the  observer. 
In  order  to  obtain  verifiable  output,  this  range  of  points  must  be  matched  to  pixels 
in  the  camera  image. 

Figure  7  shows  the  experiment  setup  showing  the  observer  (labeled  A)  the  beam 
source  (labeled  B),  the  target  (labeled  C)  and  some  point  of  interest  along  the  beam 
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(labeled  D).  We  can  obtain  the  distance  BD  using  the  Law  of  Sines  and  solving  the 
resulting  system  to  obtain 


BD  =  (AB)sin  a 
sin  d 

We  can  obtain  the  angle  b  from  the  locations  of  A,  B,  and  C,  and  the  angle  a  can 
be  obtained  directly  from  the  images  as  long  as  the  source  is  visible  in  the  image  in 
question.  The  d  is  then  simply  the  difference  n  —  a  —  b.  The  distance  AB  can  be 
determined  from  the  locations  but  laser  rangefinder  measurements  are  also  available. 

In  practice,  it  may  be  more  convenient  to  pick  the  distance  BD  and  compute  the 
angle  at  A  from  it.  Again  by  the  Law  of  Sines: 


sin  a  = 


(BD)  sin  b 
(AD) 


where  the  distance  AD  can  be  found  using  the  Law  of  Cosines: 


(25) 


AD  =  ^/(AB)2  +  (BD)2  -  2(AB)(BD)  cos  f>.  (26) 

By  a  small-angle  approximation,  each  pixel  in  the  image  represents  a  fixed  angle 
in  the  camera’s  field  of  view.  Using  this  assumption,  the  angle  a  can  be  measured 
directly  from  the  image  as  the  number  of  pixels  from  the  source  location  to  the  pixel 
of  interest.  In  this  way,  each  pixel  that  is  along  the  laser  beam  in  the  image  can  be 
related  to  a  specific  distance  from  the  source,  enabling  the  irradiance  from  that  pixel 
to  be  attributed  to  scattering  that  occurred  at  that  distance  along  the  beam. 

3.2.4  Relating  image  irradiances  to  HELEEOS  irradiances. 

HELEEOS  calculates  scattering  in  a  cylindrical  region  along  the  beam  path,  giving 
aggregate  scattering  across  the  entire  width  of  the  beam  in  a  range  of  distances 
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specified  by  the  user.  Pixel  values  in  the  image,  on  the  other  hand,  represent  the 
irradiance  produced  by  all  points  in  their  field  of  view.  Figure  9  illustrates  this. 
Typically,  a  pixel’s  field  of  view  does  not  span  the  entire  width  of  the  beam,  so  the 
irradiances  from  multiple  pixels  must  be  added  together.  A  single  irradiance  value 
can  be  computed  for  an  image  by  adding  up  pixels  covering  some  region  of  interest, 
which  can  then  be  compared  with  predictions  from  HELEEOS.  The  implementation 
of  this  process  will  be  discussed  in  more  detail  in  the  next  chapter. 


Figure  9.  HELEEOS  computes  the  total  irradiance  from  a  cylindrical  length  of  the 
beam,  while  each  pixel  of  an  image  records  the  total  irradiance  scattered  from  points 
in  its  field  of  view.  To  compare  the  two,  irradiances  from  multiple  pixels  must  be 
combined  such  that  they  contain  the  cylindrical  region  being  integrated  in  HELEEOS. 
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IV.  Analysis  Process 


Before  analyzing  the  field  images,  the  camera  had  to  be  calibrated  against  a 
known  radiation  source,  in  this  case  a  blackbody.  This  calibration  process  cleans 
up  the  images,  but  more  importantly  it  converts  the  arbitrary  units  of  the  camera 
into  physical  units  of  energy.  The  process  also  revealed  some  quirks  of  the  camera, 
and  provided  a  means  to  quantify  some  of  the  uncertainties  in  the  final  result. 

4.1  Camera  calibration 

The  camera,  a  XenICs  Xeva  409,  was  calibrated  against  a  Cl  Systems  SR80-12HT 
blackbody  at  a  variety  of  temperatures  between  200  and  500  °C,  with  the  camera 
placed  at  a  distance  of  600  cm  [16].  During  both  the  calibration  and  the  field  data 
collection,  the  camera’s  cooling  temperature  was  set  to  200  K.  The  camera  viewed 
the  blackbody  through  a  telescope  having  a  10.1  cm  aperture,  the  smaller  of  the 
two  telescopes  used  in  field  observations.  A  band-pass  filter  was  placed  in  front  of 
the  telescope,  which  has  a  pass  band  extending  from  1.49  /im  to  1.62  jim.  The 
field  of  view  of  telescope-camera  combination  was  sufficiently  small  that  the  12-inch 
blackbody  over-filled  the  detector’s  field  of  view.  Were  the  detector  not  over-filled,  a 
more  complicated  calibration  process  would  be  required. 

The  larger  of  the  two  telescopes,  a  reflecting  telescope  with  a  31.75  cm  aperture, 
needed  to  be  much  farther  away  from  its  target  in  order  to  focus,  and  adequate  space 
was  not  available  in  the  laboratory.  Therefore  the  calibration  of  the  smaller  telescope 
had  to  be  applied  to  the  larger  telescope  as  well.  This  was  possible  to  do  because 
the  same  camera  was  moved  back  and  forth  between  the  two  telescopes.  The  process 
by  which  the  calibration  results  were  applied  to  the  larger  telescope  is  discussed  in 
Section  4. 1.2.1.  Both  telescopes  and  the  camera  are  shown  in  Figure  10. 
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Figure  10.  The  telescopes  and  camera  used  for  the  field  data  collection.  The  10.10  cm 
refracting  telescope  can  be  seen  attached  to  the  side  of  the  larger  31.75  cm  reflecting 
telescope,  with  the  camera  (the  blue  cube)  attached  to  the  back  of  the  smaller  telescope. 

4.1.1  Analysis  of  the  calibration  data. 

This  section  discusses  the  process  of  analyzing  the  blackbody  images.  A  software 
application  was  developed  to  assist  with  this  analysis  and  to  apply  the  calibration 
results  to  images  taken  in  the  field.  This  program  is  documented  in  Appendix  A. 

4. 1.1.1  Spatial  calibration. 

The  first  blackbody  images  were  taken  with  a  mask  placed  over  the  blackbody. 
This  mask  provided  a  means  to  focus  the  camera,  as  well  as  to  perform  a  spatial 
calibration.  The  mask  consisted  of  a  foam-core  board  with  holes  cut  in  it  so  that  the 
blackbody  could  be  viewed  through  the  holes.  Examples  of  these  images  are  shown 
in  Figure  11.  The  dimensions  of  the  holes  were  measured  in  physical  space  with  a 
ruler,  as  well  as  in  image  space.  The  quotient  of  the  physical  extent  of  a  hole  divided 
by  the  600-cm  distance  between  the  camera  and  the  blackbody,  using  a  small-angle 
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approximation,  is  equal  to  the  angular  extent  of  the  hole.  The  ratio  of  this  angle 


to  the  hole’s  extent  in  pixels,  the  per-pixel  field  of  view  of  the  camera,  provides  a 
mapping  of  dimensions  in  image  space  to  physical  dimensions  of  the  subject  being 
imaged. 
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Figure  11.  Blackbody  images  taken  with  a  mask  over  the  blackbody,  used  for  spatial 
calibration. 


4. 1.1. 2  Selecting  blackbody  images. 

Of  the  254  blackbody  images  taken,  only  162  were  used  to  compute  the  final 
calibration.  Of  those  discarded,  a  handful  had  been  taken  with  a  mask  in  front  of 
the  blackbody  in  order  to  facilitate  spatial  calibration  as  discussed  in  the  previous 
section,  but  most  were  either  overexposed,  underexposed,  or  otherwise  flawed.  This 
section  describes  the  process  by  which  suitable  blackbody  images  were  selected. 

Histograms  of  the  blackbody  images  provided  one  means  of  assessing  the  quality 
of  a  histogram.  Since  the  scene  in  front  of  the  camera  is  mostly  a  constant  flux 
with  a  few  perturbations  added  to  it,  a  good  histogram  of  a  blackbody  looks  roughly 
Gaussian.  Before  calibration,  variations  in  the  detector  properties  over  the  surface  of 
the  focal-plane  array  (the  part  of  the  camera  that  detects  incident  photons)  widen  the 
Gaussian  distribution,  but  a  good  blackbody  image  will  retain  a  Gaussian  appearance 
both  before  and  after  calibration.  Truncation  of  the  Gaussian  shape  is  an  indication  of 
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a  large  number  of  over-exposed  or  under-exposed  pixels,  so  images  having  a  truncated 
Gaussian  distribution  were  excluded  from  the  calibration.  Examples  of  these  are 
shown  in  Figure  12. 

One  common  cause  of  overexposure  or  underexposure  is  the  use  of  an  integration 
time  that  is  too  high  or  too  low.  Integration  time  is  equivalent  to  exposure  time 
in  film  photography,  so  a  longer  integration  time  allows  more  photons  to  enter  the 
detector.  Setting  too  high  an  integration  time  will  overexpose  the  image,  while  a  low 
integration  time  can  underexpose  the  image. 

One  would  expect  the  camera  output  to  vary  linearly  with  integration  time,  but 
this  was  not  the  case,  as  seen  in  Figure  13.  The  camera’s  output  is  a  linear  function 
of  integration  time  between  about  3000  and  10000  /is,  and  a  small  region  of  linearity 
(albeit  off-set  from  the  first)  can  be  seen  between  2000  and  2500  //.$.  Between  2500 
and  3000  fis ,  however,  the  camera  outputs  varies  erratically,  with  several  points  at 
which  the  camera  outputs  zeros  throughout  the  image.  Some  of  the  large  variations  in 
camera  output  shown  in  Figure  13  occured  after  changing  the  integration  time  by  only 
1  /j,s,  and  did  so  instantly  when  the  setting  is  changed.  The  erratic  outputs  did  not 
fluctuate  substantially  with  time  (fluctuations  over  time  appeared  no  larger  than  were 
observed  for  outputs  in  the  linear  region).  Furthermore,  the  erratic  outputs  appeared 
to  be  repeatable,  such  that  when  one  observed  an  output  at  a  given  integration 
time,  changed  the  integration  time  and  then  changed  it  back,  the  output  appeared 
the  same  as  was  previously  observed  at  the  first  integration  time.  Data  were  not 
collected  to  substantiate  this  observation  about  repeatability,  but  the  abrupt  changes 
in  output  as  the  integration  time  was  changed,  combined  with  the  repeatability  of 
these  erratic  values,  suggest  that  the  behavior  was  a  quirk  of  the  camera  as  opposed 
to  the  experimental  setup. 

Images  for  which  the  output  varied  significantly  from  linearity  were  discarded  and 
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Figure  12.  Histograms  of  raw  camera  data.  Top:  Histogram  of  an  image  with  many 
underexposed  pixels.  Middle:  Histogram  of  an  image  with  many  overexposed  pixels. 
Bottom:  A  good  histogram. 
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not  used  in  the  calibration.  The  calibration  process  assumed  linearity  throughout  and 
did  not  account  for  the  offset  between  the  two  linear  regions  seen  in  Figure  13,  but 
the  error  introduced  by  this  assumption  was  captured  by  the  error  analysis  process 
described  in  Section  4.4. 


Figure  13.  Camera  response  versus  integration  time  for  a  300  C  blackbody.  Note  the 
erratic  behavior  between  2000  and  3000  /is.  This  image  shows  data  for  the  pixel  at 
(91,118),  but  other  pixels  behave  similarly. 


A  number  of  blackbody  images  showed  artifacts  of  unknown  origin,  most  com¬ 
monly  horizontal  bands  across  the  image.  These  images  were  nearly  always  outliers 
when  the  camera  output  was  plotted  against  integration  time  and  so  were  discarded. 
Examples  of  these  are  shown  in  Figure  14.  Some  of  the  blackbody  images,  partic¬ 
ularly  those  for  the  higher  temperatures,  contained  fluctuations  caused  by  hot  air 
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convecting  upward  in  front  of  the  blackbody.  These  images  were  used,  and  the  effects 
of  convection  mitigated  by  time-averaging  the  data  before  performing  the  regression. 
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Figure  14.  Blackbody  images  showing  banding  artifacts. 


4. 1.1. 3  Computing  the  calibration. 

The  calibration  was  computed  by  performing  multiple  linear  regression,  where  the 
pixel  values  from  the  camera  served  as  the  independent  variable  and  the  per-pixel  en¬ 
ergies  served  as  the  dependent  variable.  The  per-pixel  energies  were  calculated  using 
the  Planck  radiation  law  for  the  respective  blackbody  temperatures  and  integration 
times.  By  using  linear  regression,  gain  and  offset  arrays  could  be  computed  that 
provided  a  best  fit  for  all  the  blackbody  images. 

Planck  integration.  Since  the  camera  has  a  fairly  broad  wavelength  range, 
a  bandpass  filter  was  used  to  reduce  the  bandwidth  detected.  To  obtain  the  radiant 
power  output  of  the  blackbody  at  each  temperature,  the  spectral  emissivity  given 
by  the  Planck  theory  (Equations  12  and  14)  was  integrated  over  the  bandwidth  of 
the  filter.  The  Planck  emissivity  curves  are  shown  in  Figure  15.  The  emissivities 
of  the  various  curves  in  Figure  15  span  several  orders  of  magnitude.  With  these 
temperatures,  and  the  wide  variety  of  integration  times  used,  the  calibration  spanned 
the  entire  dynamic  range  of  the  camera. 
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Figure  15.  Planck  emissivity  curves  for  the  various  blackbody  temperatures. 

To  better  capture  the  filter’s  transmittance,  the  spectral  radiance  was  multiplied 
by  the  filter’s  spectral  transmittance  at  each  wavelength.  Since  the  actual  transmit¬ 
tance  was  integrated  along  with  the  spectrum,  the  integration  was  extended  well  into 
the  tails  of  the  filter  spectrum,  as  shown  in  Figure  16.  For  the  upper  limit,  the  inte¬ 
gration  was  cut  off  at  1.67  fim,  the  point  at  which  an  InGaAs  camera’s  responsivity 
falls  off  by  50%  according  to  the  responsivity  curve  in  Figure  5.  This  is  consistent 
with  the  manufacturer’s  advertisement  that  the  camera’s  spectral  band  extends  to  1.7 
/im  [20]. 1  Since  both  1.67  and  1.7  are  in  the  tail  of  the  filter  spectrum,  the  difference 
between  1.67  and  1.7  [im  will  probably  not  change  the  calibration  results  to  a  large 
degree. 

1The  Xeva  409  camera  is  no  longer  sold  and  a  specification  sheet  was  not  available;  the  source 
cited  here  is  the  specification  sheet  for  the  currently  manufactured  XS-1. 7-320,  which  is  believed  to 
be  a  nearly  identical  product. 
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Figure  16.  The  transmittance  spectrum  of  the  filter. 

Handling  bad  pixels.  A  handful  of  pixels  produced  the  same  numerical 
value  regardless  of  the  input  energy,  and  for  these  pixels  the  regression  model  failed. 
These  were  identified  as  “bad”  pixels  and  their  coordinates  stored  in  an  array. 

Once  the  regression  was  completed,  it  could  then  be  applied  to  the  blackbody  im¬ 
ages  used  in  the  calibration.  Ideally  the  calibrated  images  should  consist  of  randomly 
varying  values  with  an  average  near  the  energy  calculated  for  the  corresponding  black- 
body  temperature  and  integration  time.  Any  pixels  deviating  far  from  the  average 
for  the  image  could  be  presumed  to  have  a  significantly  nonlinear  response  and  there¬ 
fore  be  identified  as  “bad”  pixels.  The  calibration  program  provided  two  methods  for 
handling  this.  The  first  method  was  to  automatically  identify  pixels  as  bad  if  their 
value  was  more  than  some  user-specified  threshold  from  the  image’s  mean  value;  this 
threshold  was  expressed  as  a  multiple  of  the  standard  deviation  so  that  it  could  be 
applied  uniformly  to  all  the  images  in  use.  The  second  method  (used  alone  or  in  com- 
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bination  with  the  first)  was  to  allow  the  user  to  pick  pixels  to  be  treated  as  “bad”  and 
manually  add  them  to  the  array.  In  the  end,  the  automatic  method  proved  sufficient 
and  the  manual  method  was  not  used. 

When  applying  the  calibration  to  an  image,  the  pixel  values  were  multiplied  by 
the  computed  gain  array  and  then  added  to  the  computed  offset  array.  The  “bad” 
pixels  were  then  handled  by  replacing  their  values  with  the  average  of  nearby  pixels. 
Normally  this  meant  simply  averaging  all  the  good  pixels  adjacent  to  them,  but  for 
pixels  with  no  adjacent  pixels  the  program  automatically  increased  the  radius  until 
at  least  one  good  pixel  was  available  to  compute  an  average. 

Once  a  calibration  was  computed,  it  was  tested  against  some  of  the  blackbody 
images  not  used  in  the  calibration  (those  with  many  over-  or  under-exposed  pixels)  to 
verify  that  the  calibration  produced  energies  near  those  predicted  by  Planck  theory  for 
the  respective  blackbody  temperatures.  Finally,  it  was  tested  against  images  taken  in 
the  field  to  ensure  that  they  appeared  qualitatively  correct.  These  last  tests  provided 
an  additional  opportunity  to  identify  “bad”  pixels. 

4.1.2  Calibration  results. 

The  calibration  process  revealed  some  interesting  characteristics  of  the  camera. 
The  camera  has  the  characteristic  that  its  output  varies  nonlinearly  (and  sometimes 
erratically)  as  a  function  of  integration  time.  Its  response  increases  linearly  for  a  large 
range  of  integration  times,  but  with  a  few  notable  exceptions,  particularly  between 
2500  and  3000  /js,  as  shown  in  Figure  13.  These  nonlinearities  appear  in  an  identical 
fashion  when  the  camera  integration  time  is  manipulated  up  or  down;  this  repeata¬ 
bility  suggests  that  they  are  a  quirk  of  the  camera  rather  than  a  consequence  of  the 
experimental  set-up.  Images  taken  using  integration  times  within  the  range  of  linear 
response  should  be  unaffected  by  these  problems. 
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Once  the  bad  integration  times  are  removed  from  the  data  set,  Planck  theory  can 
be  used  to  compute  the  energy  expected  for  each  blackbody  temperature,  as  shown 
in  Figure  17.  The  results  show  a  weak  linear  trend  overall,  but  the  points  associated 
with  each  blackbody  temperature  are  grouped  in  lines,  and  the  points  spread  along 
these  lines  considerably  according  to  their  respective  integration  times.  Clearly  some 
significant  factor  has  not  yet  been  accounted  for  here. 


Figure  17.  Camera  response  versus  energy  received  for  various  blackbody  temperatures 
and  integration  times.  Linear  fits  are  shown  for  each  of  the  blackbody  temperatures 
for  which  at  least  2  points  are  available;  the  positive  slope  of  each  line  is  a  result  of 
increasing  integration  time.  The  range  of  integration  times  used  is  listed  in  the  legend 
of  Figure  18. 


By  grouping  points  according  to  integration  time  and  performing  linear  regression 
on  each  group,  as  shown  in  Figure  18,  it  is  seen  that  the  offset  varies  with  integration 
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time  but  the  gain  changes  relatively  little.  Plotting  the  offset  against  the  integration 
time  reveals  a  predominantly  linear  relationship,  as  shown  in  Figure  19.  A  possible 
physical  explanation  for  this  is  that  the  focal  plane  array  is  receiving  energy  from 
some  other  source  besides  the  aperture,  and  that  this  energy  arrives  at  a  roughly 
constant  rate.  This  could  be  caused  by  ambient  radiation  entering  the  aperture  at 
extreme  angles  and  scattering  inside  the  camera  body  before  reaching  the  focal  plane 
array,  or  by  electrical  processes  producing  currents  independently  of  any  radiation 
source.  Since  the  process  is  a  linear  function  of  integration  time,  it  is  possible  to 
remove  the  effect  by  linear  regression.  The  dependence  of  gain  on  integration  time, 
also  shown  in  Figure  19,  would  be  more  difficult  to  account  for,  but  is  fortunately 
less  pronounced. 

The  plot  of  offset  versus  blackbody  temperature  in  Figure  19  reveals  no  consistent 
trend,  but  the  plot  of  gain  versus  blackbody  temperature  shows  a  noticeable  decay 
as  temperature  increases.  Because  this  effect  is  related  to  blackbody  temperature,  it 
could  be  due  to  some  unaccounted  for  characteristic  of  the  experimental  setup  and 
not  the  camera.  On  the  other  hand,  the  lower  blackbody  temperature  might  produce 
a  larger  number  of  photons  which  are  near  the  edge  of  the  camera’s  response  range 
(1.7  microns),  resulting  in  a  nonlinear  response  at  that  temperature. 

The  calibration  program  was  modified  to  use  multiple  regression  in  order  to  ac¬ 
count  for  this  effect  of  the  integration  time.  Including  the  integration  time  as  part  of 
the  regression  process  produces  a  much  better  fit,  as  shown  in  Figure  20.  Additional 
blackbody  images  have  been  included  in  the  plot  to  verify  the  stability  of  the  fit.  The 
outlier  images  in  this  plot  are  extreme  cases  of  erratic  behavior  on  the  part  of  the 
camera  or  other  unusual  factors.  Each  outlier  image  contains  readily  visible  defects, 
such  as  a  horizontal  line  across  the  center  of  the  image,  or  a  histogram  showing  a 
truncated  distribution  such  as  those  shown  in  Figure  12.  Therefore,  the  calibration 
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Figure  18.  The  same  data  set  as  seen  in  Figure  17,  with  linear  fits  conducted  on  groups 
of  points  having  the  same  integration  time.  Most  of  the  lines  are  parallel,  suggesting 
that  the  offset  varies  as  a  function  of  integration  time  but  the  gain  varies  very  little. 
The  exception  to  this  is  an  integration  time  for  which  only  two  close-together  points 
were  available,  resulting  the  line  running  across  the  bottom  of  the  plot. 
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Figure  19.  Gain  and  Offset  versus  integration  time  and  temperature.  The  integration 
time  relationship  is  largely  linear  beyond  3000  /is  (a  larger  number  of  data  points  below 
2500  {is  might  reveal  a  separate  linear  domain  in  that  region),  while  the  gain  shows 
only  weak  dependence  on  integration  time.  The  offset  shows  no  consistent  temperature 
dependence,  but  the  gain  shows  a  nonlinear  temperature  dependence. 
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is  likely  to  be  reasonably  stable  when  applied  to  the  field  images. 


Blackbody  energy  received  per  pixel  (J) 


Figure  20.  Result  of  including  integration  time  in  the  regression  process. 

Figure  21  shows  the  calibration  result  across  the  focal  plane  array  in  the  form  of 
four  image  plots.  The  most  interesting  of  these  is  the  gain  array,  which  shows  a  smooth 
decrease  of  about  10%  from  outside  to  inside.  Removing  a  small- angle  approximation 
from  the  calibration  software  did  not  fix  this.  According  to  Greg  Smith,  however,  this 
particular  blackbody  exhibits  a  temperature  gradient  of  about  1.5°C,  with  the  highest 
temperatures  typically  toward  the  center  of  the  blackbody  surface.  He  said  that  the 
temperature  spread  was  1.5°C  regardless  of  the  blackbody  temperature;  it  did  not 
increase  or  decrease  substantially  when  the  blackbody  was  heated.  This  temperature 
gradient  is  a  likely  cause  of  the  main  variation  in  the  gain  array.  Mr.  Smith  allowed 
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the  blackbody  to  sit  15-20  minutes  after  the  thermostat  stabilized  before  making 
this  measurement  [17].  Many  of  the  images  used  in  the  calibration  were  taken  much 
sooner  after  the  thermostat  stabilized,  so  the  temperature  gradient  contributing  to  the 
spread  in  the  gain  array  may  be  larger  than  1.5°C.  Although  the  exact  temperature 
gradient  was  not  measured  when  the  images  were  taken,  the  errors  contributed  to  the 
calibration  by  this  gradient  can  be  captured  by  taking  the  standard  deviation  of  the 
gain  array,  since  the  temperature  gradient  appears  to  have  been  the  dominant  source 
of  variation. 

Both  the  gain  array  and  the  integration  time  coefficient  array  show  some  interest¬ 
ing  streak  features  near  the  edges  of  the  images.  The  cause  of  these  is  unknown,  but 
they  appear  not  to  be  representative  of  the  camera’s  response  in  the  field,  since  the 
same  streaking  artifacts  appear  in  some  of  the  field  images  after  calibration.  These 
features  have  magnitudes  on  the  same  order  of  the  dark  region  attributed  to  the  tem¬ 
perature  gradient,  so  the  standard  deviation  of  the  gain  array  will  also  include  errors 
contributed  by  this  effect. 

Of  the  three  arrays  used  to  apply  the  calibration,  the  offset  array  is  the  one  that 
contributes  most  to  the  appearance  of  the  output.  The  vertical  lines  seen  in  the  offset 
array  are  visible  in  every  image  produced  by  the  camera,  and  they  go  away  when  the 
calibration  is  applied. 

The  variance  array  is  computed  from  the  residual  terms  associated  with  the  linear 
regression.  It  consists  of  variances 
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where  N  is  the  number  of  images  used  in  the  calibration,  m  is  the  number  of  fit 
parameters  (3),  yimg  is  the  value  of  that  pixel  in  each  blackbody  image,  and  y  is  the 
value  predicted  by  the  calibration  for  that  blackbody  image  (from  its  temperature 
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Figure  21.  Image  plots  of  the  calibration  result.  Magenta-colored  pixels  are  considered 
“bad.”  The  gain  array  (upper  left)  characterizes  the  way  the  camera  output  varies  as 
the  incoming  radiation  is  increased.  The  dark  region  in  the  middle  of  the  gain  array 
is  probably  caused  by  a  temperature  gradient  over  the  surface  of  the  blackbody  (the 
blackbody  tends  to  be  hottest  in  the  middle).  The  offset  array  (upper  right)  represents 
the  camera’s  output  in  the  absence  of  incident  radiation;  the  vertical  lines  in  the  offset 
array  are  responsible  for  removing  similar  features  from  images  taken  with  the  camera. 
The  array  of  integration  time  coefficients  (lower  left)  is  a  correction  factor  which  had  to 
be  added  to  the  offset  array,  probably  due  to  dark  currents  or  stray  radiation  entering 
the  system.  The  variance  array  (lower  right)  is  the  sum  of  the  residuals  for  all  of  the 
images;  it  gives  an  indication  of  the  errors  associated  with  the  linear  fit. 
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and  integration  time).  This  is  a  standard  way  of  quantifying  uncertainties  for  linear 
regression  [5].  Since  the  regression  was  performed  separately  on  each  pixel,  this 
produces  an  independent  value  for  each  pixel. 

4. 1.2.1  Applying  the  calibration  data  to  a  different  optical  setup. 

Since  the  larger  telescope  could  not  be  used  to  image  the  blackbody  directly,  the 
calibration  data  from  the  small  telescope  had  to  be  applied  to  images  taken  with  the 
large  telescope  as  well.  In  order  to  do  this,  the  following  assumptions  were  made: 

•  Losses  due  to  the  optics  (from  absorption,  scattering,  transmission  through 
mirrors,  reflection  off  lenses,  etc.)  are  equal  in  both  systems. 

•  The  relative  obscuration  (the  fraction  of  the  field  of  view  that  was  blocked  by 
optical  components  such  as  mirrors)  of  the  smaller  telescope  is  zero. 

•  The  secondary  mirror  of  the  large  telescope  was  the  only  source  of  obscuration 
in  that  system  (i.e.  the  spider  and  other  parts  contributed  negligibly  to  the 
relative  obscuration). 

Using  these  assumptions,  the  calibration  could  be  applied  directly  to  the  large 
telescope,  and  the  only  change  required  was  to  use  the  larger  telescope’s  aperture 
diameter  (31.75  cm)  in  converting  from  energy  to  irradiance,  and  to  account  for  the 
relative  obscuration  by  reducing  the  aperture  area  accordingly. 

4.2  Background  subtraction 

In  the  field,  the  camera  receives  light  from  a  variety  of  sources  in  the  scene, 
and  the  laser’s  contribution  to  the  final  image  must  somehow  be  separated  from 
everything  else.  This  is  done  by  quantifying  the  background  radiation  in  the  scene 
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and  subtracting  it  from  the  calibrated  image  to  find  the  radiation  scattered  by  the 
laser. 

Many  of  the  scenes  consist  of  animations  in  which  the  laser  turns  on  or  off.  In 
this  case  the  background  can  be  computed  as  the  average  of  some  of  the  frames  in 
which  the  laser  is  turned  off.  Subtracting  this  average  from  a  frame  in  which  the  laser 
is  on  yields  an  image  in  which  most  pixels  have  values  near  zero,  except  for  those 
containing  laser  beam  scatter. 

Unfortunately,  the  background  fluctuates  slightly  with  time,  due  to  variances  in 
solar  illumination  (due  to  passing  clouds)  or  other  factors.  In  order  to  eliminate 
this  bias,  a  histogram  is  computed  from  the  image  being  studied.  Since  most  of  the 
pixels  are  background  pixels,  the  maximum  of  the  histogram  indicates  the  overall 
shift  in  background  compared  with  the  previously-computed  average.  Subtracting 
this  maximum  value  from  the  entire  image  will  remove  most  of  the  temporal  variance 
in  the  background. 

Much  of  the  background  contains  noise  which  cannot  be  eliminated.  After  shifting 
the  background,  however,  the  noise  has  an  average  value  of  zero,  so  that  positive  noise 
will,  on  average,  cancel  out  negative  noise  when  a  sum  is  computed  over  a  region  of 
the  image. 

4.3  Total  ir radiance  for  an  image 

The  calibration  process  produces  an  irradiance  for  each  pixel  of  an  image,  but 
HELEEOS  produces  irradiances  for  a  particular  span  of  the  beam  length,  which  is 
specified  as  a  distance  range  under  “Observer  Field  of  View,”  as  shown  in  Figure  22. 
In  order  to  compare  HELEEOS  outputs  to  the  measured  irradiances  in  a  meaningful 
way,  this  distance  range  must  be  mapped  to  a  set  of  pixels  in  each  image.  Once 
the  appropriate  pixels  are  selected,  the  total  irradiance  for  the  region  can  be  found 
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by  adding  together  the  irradiances  for  all  the  pixels  in  the  region.  Ideally  this  total 
should  be  close  to  the  irradiance  produced  by  HELEEOS  as  long  as  the  distance  range 
and  other  parameters  are  consistent. 


Figure  22.  Observer  optics  settings  in  HELEEOS.  The  analysis  software  maps  the 
distance  ranges  under  “Linear  Focusing”  (here  identified  by  a  red  box)  to  a  region  of 
pixels  in  the  image. 


For  each  image,  a  range  of  distances  along  the  laser  beam  is  selected.  The  dis¬ 
tances,  defined  in  meters  from  the  laser  source,  are  chosen  such  that  the  entire  range 
is  visible  in  the  image.  This  distance  range  is  converted  to  angles  using  Equation  25, 
and  these  angles  are  converted  to  distances  in  pixels  using  the  per-pixel  field  of  view 
computed  in  Section  4. 1.1.1. 
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The  selection  region  is  defined  in  terms  of  the  location  of  the  laser  source  in  the 
image,  the  coordinates  of  which  are  found  manually  for  each  image.  The  selection 
region  is  defined  relative  to  the  source,  as  shown  in  Figure  23.  The  parameters  are 
chosen  in  such  a  way  as  to  contain  the  entire  beam  within  the  desired  distance  range. 
The  region  is  given  in  terms  of  an  apparent  beam  direction,  0,  which  defines  the 
direction  of  the  beam  in  the  image,  clockwise  from  the  horizontal  direction.  The 
beam’s  width  in  the  image  is  defined  in  terms  of  a  parameter  A0,  which  specifies  a 
range  of  angles  on  either  side  of  the  beam’s  center. 


Figure  23.  Parameters  defining  the  selection  region.  All  angles  and  distances  are 
defined  relative  to  the  location  of  the  beam  source,  and  were  entered  manually.  The 
angle  0  defines  the  apparent  direction  of  the  beam,  clockwise  from  the  horizontal  axis. 
The  selection  region  includes  all  pixels  whose  direction  from  the  beam  source  is  between 
6  — AO  and  6  + AO,  and  whose  distance  from  the  source  is  FOV  Start  and  FOV  Stop.  FOV 
Start  and  FOV  Stop  are  typically  entered  in  meters  (representing  distances  along  the 
beam  from  the  source),  and  are  converted  to  angles  and  pixel  ranges  in  the  observer’s 
field  of  view  using  Equation  25. 


The  total  irradiance  for  the  image  is  found  by  summing  the  irradiances  computed 
for  each  pixel  within  the  selected  region.  This  summation  is  computed  for  each  frame 
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within  a  range  of  interest.  The  frame  range  is  selected  such  that  the  laser(s)  are 
turned  on  and  in  a  steady  state  (i.e.  they  are  not  in  the  middle  of  turning  on  or 
off).  Sometimes  several  frame  ranges  were  selected  out  of  the  same  data  set.  This 
was  typically  the  case  when  the  lasers  were  turned  on  or  off  one  by  one,  so  that  part 
of  the  data  set  showed  one  laser  turned  on,  another  showed  two  lasers,  and  perhaps 
another  section  showed  three  lasers.  For  each  frame  range,  the  irradiance  reported 
as  the  final  result  was  the  average  of  the  totals  for  all  the  frame,  with  the  standard 
deviation  of  those  totals  being  used  as  an  error  term. 

4.4  Error  analysis 

Known  and  quantifiable  errors  in  the  analysis  include:  Random  variances  in  the 
calibration  data,  nonuniformity  of  the  blackbody  temperature  over  its  surface,  non¬ 
zero  background  values  that  were  not  removed  by  background  subtraction,  and  time- 
varying  fluctuations  in  the  irradiance  from  one  frame  to  another. 

Random  differences  between  blackbody  images  are  computed  using  the  expression 
for  a  fa,  given  in  Equation  27.  afit  includes  errors  due  to  the  blackbody  temperature 
being  higher  or  lower  than  it  was  assumed  to  be,  changes  in  the  energy  entering  the 
camera  due  to  variations  in  room  lighting,  and  changes  in  the  state  of  the  camera 
over  the  course  of  the  calibration.  It  also  includes  errors  introduced  by  the  assump¬ 
tion  that  the  camera’s  response  is  a  linear  function  of  incident  energy,  which  is  of 
course  a  simplification.  The  assumption  that  photon  count  and  energy  are  equivalent 
introduced  a  small  amount  of  error  which  is  captured  in  a f  lt .  Errors  in  measuring  the 
distance  between  the  aperture  and  the  blackbody  are  also  effectively  included  in  cr/*t, 
since  the  calibration  was  conducted  over  the  course  of  three  days  so  this  distance  was 
measured  several  times.  Even  so,  the  600  cm  distance  is  considered  to  be  accurate  to 
within  1  cm;  at  0.2%,  this  is  smaller  than  most  of  the  other  errors  being  considered. 
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As  mentioned  previously,  the  dark  area  in  the  middle  of  the  gain  array  was  caused 
by  a  temperature  gradient  from  the  center  to  the  edges  of  the  blackbody.  In  addition, 
there  are  horizontal  and  vertical  streaks  near  the  edges  of  the  gain  array.  The  cause 
of  these  streaks  is  unknown,  but  they  may  well  be  caused  by  errors  in  the  calibration 
process  (such  as  an  additional  region  of  temperature  variation  on  the  blackbody 
surface),  rather  than  an  actual  variation  in  gain  from  one  pixel  to  another.  In  any 
case,  these  streaking  features  introduce  gain  variations  on  the  order  of  those  caused 
by  the  temperature  fluctuations,  as  evidenced  by  the  fact  that  both  appear  blue  in 
Figure  21.  Therefore,  the  standard  deviation  og  of  the  gain  array  can  be  considered 
to  capture  the  impact  of  the  blackbody  temperature  variations  and  other  errors  upon 
the  gain  array. 

The  temperature  displayed  on  the  blackbody  thermostat  comes  from  a  thermo¬ 
couple  located  behind  the  center  of  the  blackbody  surface.  The  temperature  from 
this  thermocouple  is  certified  by  the  manufacturer  to  be  accurate  to  within  0.5°C. 
The  manufacturer’s  acceptance  test  report  lists  actual  errors  measured  in  testing  the 
blackbody,  all  of  which  were  considerably  smaller  than  0.5°C  (they  ranged  from  0.01 
to  0.07°C)  [7].  Deviations  of  the  blackbody  temperature  would  bias  the  result,  but 
since  the  0.5°C  is  less  than  the  1.5°C  variation  in  temperature  over  the  blackbody 
surface  (which  has  already  been  quantified) ,  the  uncertainty  of  the  thermocouple  was 
not  included  in  the  error  analysis. 

After  performing  background  subtraction,  most  pixels  of  the  image  should  have 
a  value  near  zero,  with  only  a  handful  containing  the  signal.  Therefore  the  standard 
deviation  cq,s  of  all  the  pixel  values  in  the  image  (after  background  subtraction)  should 
be  representative  of  the  errors  contributed  by  background  subtraction. 

Following  the  standard  procedure  of  summing  errors  in  quadrature  to  obtain  the 
uncertainty  for  a  summed  value,  the  total  uncertainty  in  the  energy  Q  received  at  a 
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pixel  is  given  by 


^Qpx 


(28) 


where  Dpx  is  the  value  in  digital  units  produced  by  the  camera  at  the  pixel  of  interest, 
minus  the  value  of  the  (un-calibrated)  background  at  that  pixel. 

The  total  energy  received  is  found  by  summing  Q  over  a  region  of  interest.  The 
uncertainties  a q  are  not  correlated,  so  they  are  summed  in  quadrature  to  obtain 
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When  the  irradiance  is  computed,  this  energy  is  multiplied  by  a  conversion  factor 
to  get  irradiance  as 


E  = 


Spx  (Q  px) 

'4-tint(i-o 


(30) 


where  d  is  the  aperture  diameter,  tvnt  is  the  integration  time,  and  C  is  the  relative 
obscuration.  The  aperture  diameter  d\  is  known  to  1  mm,  as  is  the  diameter  of  the 
obscuring  secondary  mirror  d‘> ,  which  is  responsible  for  C  in  the  larger  telescope. 
Thus  di  and  C  contribute  uncertainties  as  follows: 


dE  _  8  Zpx  (Qpx) 
dd '°dl  “  (C-l)<%nadt 
9E  _  4  Zpx  (Qpx) 
dCac  ~  (C  -  1  )2d?7T<Tc 


where  and  a c  are  the  uncertainties  in  d\  and  C,  respectively. 

C  is  equal  to  the  ratio  of  the  area  of  the  obscuring  optic  to  the  area  of  the  whole 

d2 

aperture,  that  is,  C  =  Thus  the  uncertainty  in  the  relative  obscuration  C  is 
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This  leads  to  an  expression  for  ge , 
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Variations  in  the  signal  level  over  time  may  be  caused  by  fluctuations  in  the 
background  radiation,  random  noise  due  to  the  camera  electronics,  changes  in  laser 
output,  and  variations  in  atmospheric  factors  such  as  aerosol  concentrations.  Since 
the  camera  captures  video,  the  irradiance  is  computed  separately  for  each  frame  in  a 
series  of  interest,  and  the  final  result  is  reported  as  an  average  of  these.  The  standard 
deviation  at  of  these  irradiances,  then,  provides  an  estimate  of  the  time-varying  error 
terms.  To  compute  the  uncertainty  of  the  time- averaged  values,  the  existing  ge  must 
be  broken  into  its  time-dependent  and  time-independent  parts. 

For  this  purpose,  time-independent  error  terms  will  still  be  considered  time- 
independent  if  they  are  multiplied  by  a  time-dependent  coefficient.  For  instance, 


8  Yjpx  (  Qpx ) 

-  l)d\n°dl 

is  considered  time-independent  because  a ^  is  time-independent,  even  though  Qpx  may 
vary  with  time.  Although  individual  pixels  may  fluctuate  considerably  over  time,  the 
frame  ranges  have  been  selected  such  that  the  sum  J2px  (Qpx)  should  not  fluctuate 
too  drastically.  The  only  time-dependent  error  term  which  must  be  considered  is  cp,s, 
which  can  vary  as  the  background  fluctuates.  All  the  other  error  terms  are  the  result 
of  measurements  made  prior  to  the  collection. 

When  the  images  are  averaged,  they  are  summed  and  divided  by  the  number  of 
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frames.  Since  the  background  is  simply  subtracted  from  the  image,  it  is  effectively 
summed  as  well  and  divided  by  the  same  factor.  The  fluctuating  background  uncer¬ 
tainties  do  not  correlate  with  one  another,  so  they  are  summed  in  quadrature  and 
then  divided,  so 


9avg 


S frames  ^b9  f 


rame 


N 


where  N  is  the  number  of  frames  being  averaged. 

Separating  from  the  other  terms  yields  the  time  averaged  cte 
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4.5  Computing  irradiances  in  HELEEOS 

For  each  image,  the  wavelengths  and  power  levels  of  each  laser  were  recorded 
using  the  calibration  software,  so  that  they  would  be  stored  alongside  the  rest  of  the 
data  associated  with  that  image.  A  script  was  written  which  ran  HELEEOS  with  the 
parameters  for  each  image  (and  each  laser  present  in  that  image).  When  multiple 
lasers  were  present  in  an  image,  HELEEOS  was  run  once  for  each,  and  the  total 
reported.  The  results  of  these  HELEEOS  outputs  will  be  reported  in  the  following 
section  and  compared  with  the  irradiances  measured  from  the  corresponding  images. 
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V.  Field  Measurements 


Field  measurements  were  taken  of  a  laser  test  conducted  at  Wright-Patterson 
Air  Force  Base  Area  B  in  July  2009.  The  test  was  conducted  by  AFRL  over  three 
days  with  the  primary  purpose  of  testing  a  new  technique  for  measuring  atmospheric 
turbulence.  The  three  lasers  used  in  the  test  operated  at  fixed  wavelengths  1.533, 
1.589,  and  1.595  /im,  respectively,  and  their  power  levels  were  adjusted  from  0.6  to 
5.7  W  during  the  test.  The  two  days  of  test  data  provided  results  for  two  different 
off-axis  angles,  nighttime  and  daytime  conditions,  as  well  as  a  significant  amount  of 
variation  in  laser  power.  Each  data  set  was  processed  as  described  in  the  previous 
chapter,  producing  an  irradiance  value  for  each  image.  These  irradiances  were  then 
compared  with  HELEEOS  outputs  for  the  same  geometry  and  laser  properties. 

5.1  Experiment  setup 

Lasers.  The  lasers  used  were  IPG  Photonics  erbium  fiber  lasers,  models 
ELR-5-1595.3,  ELR-5-1532.6,  and  ELR-5-1588.7.  The  spec  sheets  for  each  laser  show 
their  bandwidths  to  be  ±0.1  nm,  small  enough  to  consider  them  monochromatic  for 
this  study.  When  operating  from  an  outdoor  location,  the  beams  were  pointed  out 
the  door  of  a  trailer  as  shown  in  Figure  24.  The  lasers  were  reported  to  have  an  exit 
aperture  of  4  mm  and  a  spot  size  of  1.6  m  at  the  target,  so  these  values  were  entered 
into  HELEEOS  as  shown  in  Figure  25 

Experiment  geometry.  The  15  July  2009  test  was  conducted  with  a  source 
laser  at  the  west  end  of  the  abandoned  runway  at  Patterson  Field,  with  the  target  on 
the  ninth  floor  of  Building  622  [8].  Building  622,  including  the  window  which  served 
as  the  target,  is  shown  in  Figure  26.  The  total  distance  from  source  to  target  was 
approximately  3050  m.  The  observer  was  located  on  a  concrete  pad  on  the  west  side 
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Figure  24.  The  trailer  from  which  the  lasers  operated  when  located  out  doors.  The 
apertures  for  each  laser  are  at  the  front  of  the  black  cylindrical  object  pointing  out  of 
the  trailer. 

of  the  parking  lot  at  Building  622,  approximately  60  m  from  the  target  and  3015  m 
from  the  source.  The  distance  from  observer  to  target  was  measured  using  a  laser 
rangefinder;  all  other  distances  were  measured  using  Google  Earth.  This  setup  is 
shown  in  Figure  27. 

The  source  was  located  at  an  elevation  of  240  m,  while  the  target  was  at  313 
m.  The  observer  was  located  at  an  elevation  of  292  m.  The  observer  elevation  was 
computed  relative  to  the  building  elevation  by  analyzing  Figure  26.  The  elevation 
and  altitude  differences  altogether  place  the  observation  location  approximately  0.9 
degrees  off-axis. 

For  the  17  July  2009  test,  the  source  laser  was  located  near  the  east  end  of  the 
abandoned  runway  at  Patterson  Field,  with  the  target  2160  meters  away  at  the  west 
end  of  the  runway.  Off-axis  scattering  was  observed  from  a  location  roughly  46  meters 
from  the  target  in  a  direction  perpendicular  to  the  beam  path.  Thus  the  observer 
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Figure  25.  Optics  settings  for  the  lasers.  Non-default  values  are  marked  in  red. 
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Figure  26.  Building  622  at  Wright-Patterson  Air  Force  Base.  The  target  was  in  the 
left-most  of  the  two  visible  windows. 


Figure  27.  Experiment  geometry  for  the  15  July  2009  test.  For  this  test,  the  observer 
was  located  almost  directly  below  the  beam.  The  300  m  field  of  view  is  representative 
of  the  31.75  cm  reflecting  telescope,  which  was  used  for  most  of  the  observations  on  15 
July. 
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location  is  around  1.22  degrees  off-axis.  This  setup  is  shown  in  Figure  28.  The 
relevant  distances  were  measured  using  a  laser  rangefinder.  Observer,  source,  and 
target  were  all  located  at  roughly  the  same  elevation  of  240  m. 


Figure  28.  Experiment  geometry  for  the  17  July  2009  test.  The  source  is  denoted  by  a 
laser  symbol,  the  target  by  a  bulls-eye,  and  the  observer  by  an  eyeball.  The  observer’s 
field  of  view  is  identified  by  a  blue  line  along  the  laser  beam  path.  The  400  m  field  of 
view  is  representative  of  the  10.1  cm  refracting  telescope,  which  was  used  for  most  of 
the  observations  on  17  July. 

During  the  course  of  the  experiment  the  lasers  were  sometimes  walked  off  from 
the  target,  but  the  precise  direction  of  the  walk-off  was  not  determined.  Therefore, 
in  computing  the  observer’s  off-axis  angle  it  is  assumed  that  the  lasers  are  always 
pointed  at  the  intended  target.  This  off-axis  angle  is  used  both  in  determining  the 
observation  path  and  as  an  input  parameter  to  HELEEOS. 

Environmental  conditions.  Both  days  of  testing  were  marked  by  fair  weather, 
with  temperatures  around  20° C  and  scattered  to  broken  cloud  cover.  National  Weather 
Service  records  show  the  visibilities  for  both  days  to  be  10  statute  miles,  but  their 
measurements  were  limited  to  the  farthest  visible  object.  Those  who  were  present 
at  the  test  reported  that  the  visibilities  were  extremely  high,  far  above  the  10  miles 
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reported. 

The  meteorological  settings  entered  into  HELEEOS  were  as  follows:  The  ap¬ 
proximate  coordinates  of  Wright-Patterson  Air  Force  Base  were  selected  on  the  map 
and  the  DAYTON/WRIGHT-PATTE  (sic)  ExPERT  site  was  selected.  The  “En¬ 
able  Advanced  Atmosphere  Settings”  checkbox  was  checked,  as  was  the  “Enable  User 
Specified  Ground  Visibility,”  checkbox,  with  the  Aerosol  Visibility  setting  set  to  100 
km.  With  these  settings,  HELEEOS  was  set  to  the  ExPERT  database  for  all  weather 
parameters,  except  for  the  visibility  which  was  set  to  a  relatively  high  value  to  reflect 
the  fact  that  the  visibility  was  reported  to  be  high  by  those  present  at  the  test.  The 
HELEEOS  Scenario  pane  is  shown  in  Figure  29  with  these  settings  entered. 


Figure  29.  HELEEOS  Scenario  page,  with  fields  populated  with  the  input  data  used 
for  HELEEOS  runs  in  the  study. 
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5.2  Image  selection 


In  order  for  a  data  collection  to  be  used,  atmospheric  scattering  must  be  visible  in 
the  images,  and  it  must  be  possible  to  determine  which  of  the  lasers  is  responsible  for 
the  scattering.  In  addition,  it  must  be  possible  to  determine  the  physical  location  of 
the  beam  from  image  coordinates.  In  order  to  do  this,  images  were  selected  in  which 
the  laser  source  was  visible  in  the  picture  and  in  which  the  camera  was  stationary.1 
Finally,  a  background  image  must  be  available.  Ideally  this  can  be  created  from  an 
average  of  frames  in  the  same  data  set.  For  instance,  if  the  movie  shows  the  lasers 
being  turned  on  or  off  it  will  often  include  a  substantial  number  of  frames  in  which  no 
lasers  are  turned  on.  Otherwise,  a  background  image  can  be  constructed  from  another 
image  of  the  same  scene.  When  the  background  was  constructed  from  another  image, 
the  result  was  generally  inferior,  but  images  processed  in  this  way  were  kept  since 
they  did  not  consistently  appear  as  outliers  in  the  final  data  set. 

The  quality  of  the  images  used  varied  significantly.  Most  notable  was  the  varying 
effectiveness  of  the  background  subtraction  technique.  The  best  results  occurred 
when  the  background  could  be  computed  from  different  frames  of  the  same  hie,  while 
subtracting  a  background  computed  from  a  different  hie  yielded  far  poorer  results. 
Examples  of  the  background  subtraction  results  are  shown  in  Table  1. 

The  left-hand  column  of  Table  1  shows  a  highly  successful  background  subtrac¬ 
tion.  The  resulting  image  contains  near-zero  values  almost  everywhere  except  where 
the  laser  beam  can  be  seen.  The  center  column  shows  a  moderately  successful  back¬ 
ground  subtraction.  In  this  image,  background  features  remain  visible  in  the  image 
after  background  subtraction  has  been  performed,  but  the  laser  beam  pixels  generally 
contain  much  larger  magnitudes  than  do  background  pixels. 

1It  should  be  possible  to  handle  the  case  of  a  panning  camera  using  motion  tracking  software,  or 
with  other  equipment  it  might  have  been  possible  to  account  for  panning  by  measuring  the  camera 
direction  in  the  field. 
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Source  file 


Original 


Background 


Background- 

subtracted 


Table  1.  Examples  of  background  subtraction. 
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The  right-hand  column  of  Table  1  shows  a  relatively  unsuccessful  background 
subtraction.  Many  background  features  have  magnitudes  at  or  above  those  of  the 
laser  beam  pixels,  even  after  background  subtraction.  Still,  the  irradiances  computed 
for  these  images  did  not  tend  to  be  outliers  when  compared  with  those  computed 
from  other  images.  As  a  result,  less  successful  cases  of  background  subtraction  were 
included  in  the  final  data  set. 

5.3  Results 

In  total,  the  calibration  was  applied  to  56  images  taken  at  the  July  2009  laser 
test.  A  detailed  table  of  the  results  can  be  seen  in  Appendix  C.  For  each  image,  an 
irradiance  value  was  computed  using  HELEEOS  and  compared  with  that  measured 
from  the  image.  In  a  slight  majority  of  cases,  HELEEOS  calculated  irradiances  1-2 
orders  of  magnitude  higher  than  that  measured  in  the  field  for  the  same  scenario,  as 
shown  in  Figure  30. 

Each  image  showed  scattering  from  either  one,  two,  or  all  three  lasers.  When 
multiple  laser  beams  were  visible  simultaneously,  the  sum  of  their  power  settings  was 
shown  as  the  “Laser  power”  on  the  x-axis  of  Figure  30.  In  a  majority  of  cases,  the 
lasers  each  operated  at  0.6  W,  accounting  for  the  clusters  of  data  at  0.6,  1.2,  and  1.8 
W.  The  three  points  plotted  at  around  17  W  consisted  of  all  three  lasers  operating 
at  5. 5-5. 7  W. 

Figure  30  shows  a  weak  downward  trend  in  irradiance  ratio  as  the  laser  power 
increases.  There  is  not  enough  data  at  the  higher  power  levels  to  ascertain  whether 
this  trend  is  representative  of  a  real  phenomenon  or  whether  it  is  simply  a  side-effect 
of  the  small  size  of  the  data  set. 

Another  view  of  this  data  set  can  be  seen  in  Figure  31.  This  plot  shows  a  weak 
correlation  between  measured  irradiance  and  HELEEOS  output.  It  is  possible  that 
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Figure  30.  Ratios  of  HELEEOS-derived  irradiances  to  the  measured  irradiances,  as  a 
function  of  laser  power.  HELEEOS  outputs  are  always  within  two  orders  of  magnitude 
of  the  measured  irradiances,  and  frequently  within  one  order  of  magnitude. 
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this  correlation  would  appear  stronger  if  more  images  were  available  at  the  higher 
laser  powers,  which  would  help  show  whether  some  of  the  poorly  correlated  points  in 
this  data  set  are  anomalous  or  whether  they  are  representative  of  typical  variations. 


HELEEOS  irradiance  (W/m2  ) 


Figure  31.  Comparison  of  HELEEOS-derived  irradiances  with  measured  irradiances, 
showing  a  weak  correlation. 

In  nearly  all  of  the  field  images,  it  is  apparent  that  irradiance  decays  along  the 
length  of  the  beam.  To  determine  whether  HELEEOS  outputs  agreed  with  this  trend, 
the  the  FOV  Start  and  FOV  Stop  distances  of  the  selection  region  were  converted  to 
angles  in  the  observer’s  field  of  view,  and  the  angle  range  broken  into  increments  of 
equal  sizes,  as  shown  in  Figure  32.  The  start  and  stop  angle  of  each  increment  was 
converted  back  into  a  distance  range,  which  was  passed  to  HELEEOS  as  the  FOV 
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Starting  Dist  and  FOV  Ending  Dist  parameters.  By  running  HELEEOS  once  for  each 
angle  increment,  the  trend  of  observer  irradiance  as  a  function  of  viewing  angle  could 
be  determined.  Similarly,  the  irradiances  in  the  images  were  summed  over  each  angle 
increment,  producing  measured  irradiance  values  at  intervals  over  the  same  range  as 
was  measured  for  HELEEOS. 


Figure  32.  Angle  increments  within  the  selected  field  of  view,  used  for  computing 
irradiances  as  a  function  of  viewing  angle  for  Figure  33. 

Results  of  this  process  for  two  images  can  be  seen  in  Figure  33.  For  both  im¬ 
ages,  a  downward  trend  in  irradiance  can  be  seen  in  both  the  measurements  and  the 
HELEEOS  outputs,  although  the  overall  magnitudes  of  the  irradiances  do  not  agree 
very  well,  nor  do  the  rates  at  which  the  irradiance  falls  of  with  angle. 

The  beam  could  in  some  cases  be  detected  at  viewing  angles  as  large  as  20°,  but 
at  more  extreme  angles  it  became  too  dim  for  the  camera  to  detect.  Using  the  same 
process  used  previously  to  divide  the  field  of  view  into  equal  angle  increments,  a 
succession  of  HELEEOS  runs  was  performed  using  equal  angle  increments  from  0° 
to  90°.  In  results  of  this  test,  shown  in  Figure  34,  HELEEOS  predicts  the  irradiance 
will  continue  to  fall  off  through  75°,  which  is  consistent  with  the  fact  that  the  beam 
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0006. xvi: 


Figure  33.  HELEEOS-predicted  irradiance  (solid  lines)  and  measured  irradiance 
(points)  as  a  function  of  observation  direction  for  two  images.  Although  the  over¬ 
all  magnitude  differs  between  HELEEOS  and  the  measured  irradiance,  as  does  the  rate 
at  which  the  irradiance  falls  off  with  distance,  the  downward  trend  is  present  in  both. 
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was  invisible  to  the  camera  at  that  angle. 


Figure  34.  HELEEOS-predicted  irradiance  as  a  function  of  observation  direction  for  an 
observer  located  1.22°  off-axis.  Each  point  along  the  line  represents  an  equal  span  of 
angles  in  the  observer’s  field  of  view.  An  observation  angle  of  0°  indicates  the  observer 
is  looking  directly  at  the  source;  90°  indicates  the  observer  is  looking  90°  away  from 
the  beam  source,  approximately  perpendicular  to  the  beam.  HELEEOS  outputs  NaN 
for  the  irradiance  from  75-90°. 


The  investigation  into  the  viewing  direction  begs  another  question:  How  does  the 
observer’s  location  off-axis  affect  the  HELEEOS  output?  Since  the  Mie  phase  func¬ 
tions  show  strong  forward  scattering  and  a  fairly  rapid  fall-of  from  the  0°  direction, 
one  would  expect  that  the  smallest  off-axis  angles  (the  angle  between  the  observer 
and  the  laser  beam  path)  will  produce  the  greatest  off-axis  irradiances,  with  fairly 
rapid  fall-off  from  there.  Figure  35  shows  HELEEOS-predicted  irradiances  for  a  range 
of  small  off-axis  angles.  As  expected,  the  highest  irradiances  occur  at  the  smallest 
angles,  with  a  sharp  reduction  in  irradiance  as  the  observer  is  moved  away  from  the 
beam. 
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Figure  35.  HELEEOS-predicted  irradiance  as  a  function  of  observer  off-axis  angle. 
Larger  angles  indicate  the  observer  was  further  offset  from  the  beam  path,  maintain¬ 
ing  the  same  distance  from  the  laser  source.  The  strong  dependence  of  the  output 
on  this  parameter  suggests  that  accurate  results  from  HELEEOS  require  an  accurate 
measurement  of  the  experiment  geometry. 
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The  strong  sensitivity  of  HELEEOS  outputs  to  off-axis  angle  indicates  that  the 
off-axis  angle  must  be  determined  precisely  for  the  results  to  be  accurate.  The  off- 
axis  angles  used  for  the  test  data,  however,  could  only  be  determined  to  within  0.5°, 
but  their  magnitudes  were  on  the  order  of  1°.  This,  combined  with  imprecisions  in 
the  other  input  parameters,  provides  one  possible  explanation  for  the  wide  disparities 
seen  between  HELEEOS  and  the  test  data. 

At  worst,  HELEEOS  agrees  with  the  test  data  to  within  two  orders  of  magnitude, 
although  frequently  it  does  much  better  than  this,  in  some  cases  producing  a  value 
within  a  few  percent  of  the  test  data.  In  addition,  HELEEOS  correctly  predicts  that 
irradiance  should  decay  with  respect  to  viewing  direction.  Imprecisions  in  the  input 
parameters  might  explain  many  of  the  discrepancies.  This  test  is  not  sufficient  to 
validate  HELEEOS,  but  it  shows  that  HELEEOS  has  potential  to  be  validated  with 
further  testing. 
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VI.  Conclusions 


The  tests  conducted  in  this  research  show  some  agreement  between  HELEEOS- 
predicted  off-axis  irradiances  and  those  measured  using  a  camera.  The  data  were  not 
in  close  enough  agreement  with  HELEEOS  to  serve  as  a  validation  of  HELEEOS  by 
themselves,  but  HELEEOS  showed  some  correlation  with  the  test  data  in  terms  of 
laser  power  output  and  viewing  angle. 

For  the  entire  set  of  test  data,  HELEEOS  produced  outputs  within  two  orders  of 
magnitude,  and  frequently  within  one  order  of  magnitude.  Although  these  errors  are 
large,  there  are  cases  in  which  precision  such  as  this  would  be  adequate.  For  instance, 
if  one  wished  to  determine  whether  a  certain  laser  could  be  detected  under  a  particular 
set  of  conditions,  one  could  use  HELEEOS  to  do  this.  As  long  as  HELEEOS  predicted 
irradiances  more  than  two  orders  of  magnitude  above  or  below  the  threshold  for 
detection,  one  could  make  a  fairly  certain  determination  as  to  whether  the  laser  could 
be  detected  under  those  conditions.  For  a  laser  weapon,  a  two  order  of  magnitude 
accuracy  would  in  certain  cases  be  sufficient  to  determine  whether  a  laser  was  capable 
of  destroying  a  particular  target. 

Although  the  two  order  of  magnitude  threshold  sets  an  upper  bound  on  the  uncer¬ 
tainty  of  HELEEOS  results  for  this  set  of  operating  conditions  (clear  air,  f  .5-1.6  /j,m, 
0.6-5. 7  W,  and  camera  roughly  1°  off  axis),  the  real  value  of  this  test  is  to  demonstrate 
techniques  for  collecting  and  analyzing  off-axis  scattering  data,  and  to  provide  insights 
into  the  best  path  toward  validating  HELEEOS.  The  test  showed  that  an  infrared 
camera  can  be  used  successfully  to  measure  off-axis  scattering  data,  and  techniques 
were  developed  that  can  be  used  to  relate  that  data  to  HELEEOS  inputs  and  out¬ 
puts,  both  for  purposes  of  validation  and  for  later  real-world  applications.  The  large 
discrepancies  between  HELEEOS  outputs  and  measured  irradiances  suggest  that  a 
sensitivity  study  of  the  off-axis  scattering  model  should  be  accomplished  in  order  to 
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guide  further  testing. 

The  quality  of  the  input  data  to  HELEEOS  appears  to  have  had  a  considerable 
influence  on  the  off-axis  scattering  results.  The  off-axis  angle  (the  angle  at  which 
the  observer  was  off-set  from  the  laser  beam,  measured  from  the  source  location) 
had  magnitudes  of  1.2°  on  the  first  day  of  testing  and  0.9°  on  the  second,  but  its 
uncertainty  was  on  the  order  of  0.5°.  Furthermore,  it  has  been  shown  that  HELEEOS 
off-axis  scattering  outputs  are  extremely  sensitive  to  variations  in  this  parameter 
when  it  is  so  small,  so  that  changes  of  ±0.5°  could  change  the  results  by  an  order  of 
magnitude. 

The  atmospheric  visibility  could  also  influence  the  results  by  significant  factors. 
Although  the  outputs  are  less  sensitive  to  atmospheric  visibility  than  they  are  off- 
axis  angle,  the  range  of  possible  visibilities  was  considerable.  Atmospheric  visibilities 
are  typically  reported  based  on  the  most  distant  visible  object.  Weather  reports  for 
Wright-Patterson  Air  Force  Base  never  include  visibilities  over  10.0  miles  because 
even  on  a  perfectly  clear  day  no  objects  beyond  that  distance  can  be  seen.  Those 
present  at  the  test  reported  that  the  visibility  during  the  test  was  much  higher  than 
that,  but  they  could  not  quantify  the  actual  visibility.  A  visibility  of  100  km  was 
chosen  because  it  was  consistent  with  the  reports  of  excellent  visibility,  but  this  value 
could  probably  be  set  within  a  range  of  50-200  km  and  remain  consistent  with  those 
reports.  Changes  within  this  range  can  change  the  outputs  by  up  to  an  order  of 
magnitude. 

The  background  subtraction  did  not  always  work  as  well  as  one  would  hope, 
because  for  some  of  the  images  the  closest  available  background  image  had  been  taken 
several  tens  of  minutes  before  or  after  the  actual  data.  This  certainly  affected  the 
accuracy  of  the  irradiance  values  computed  from  these  images.  The  best  backgrounds 
were  contained  within  the  same  animation  file  as  the  test  data,  so  that  the  two  were 


69 


very  close  in  time.  Future  studies  should  take  care  to  obtain  high  quality  background 
images,  as  this  can  make  a  significant  difference  in  the  quality  of  results. 

Another  possible  source  of  error  is  the  laser  power  output.  The  power  settings 
passed  to  HELEEOS  came  from  the  amount  of  electrical  power  used  by  the  lasers,  but 
in  practice  the  optical  power  produced  by  the  laser  is  probably  significantly  less  than 
this.  Placing  a  power  meter  in  the  beam  path  might  help  correct  for  this,  although 
having  one  there  continuously  would  partially  block  the  beam  and  reduce  the  amount 
of  scattered  radiation. 

Unfortunately,  the  data  available  gives  no  immediate  indication  as  to  whether  the 
discrepancies  are  due  to  errors  in  HELEEOS  or  due  to  problems  with  the  input  data. 
Nonetheless,  it  gives  some  indications  that  the  HELEEOS  off-axis  model  is  at  least 
qualitatively  accurate.  The  data  also  shows  that  HELEEOS  irradiance  varies  with 
viewing  angle  in  a  manner  similar  to  the  physical  reality.  HELEEOS  correctly  predicts 
that  irradiance  should  decay  with  viewing  angle,  and  at  extreme  angles  it  produces 
very  low  irradiances  consistent  with  the  fact  that  the  beam  was  not  observable  at 
those  angles. 

By  collecting  and  analyzing  off-axis  scattered  radiation,  it  has  been  demonstrated 
that  off-axis  scattering  can  be  quantified  for  communications  lasers,  suggesting  that 
with  further  refinement  of  the  HELEEOS  model  and  more  accurate  input  parameters, 
much  more  consistent  results  can  be  obtained. 

6.1  Future  work 

The  existing  data  provides  ample  opportunities  for  additional  work.  Many  of  the 
images  available  were  not  used  because  of  flaws  which  could  be  overcome  with  im¬ 
provements  to  the  calibration  software.  HELEEOS  runs  compared  with  these  data 
were  done  using  climatological  weather  and  mostly  default  weather  settings,  but 
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studies  could  be  done  using  actual  weather  conditions  as  well.  Sensitivity  studies 
of  HELEEOS  would  help  determine  which  input  parameters,  meteorological  or  oth¬ 
erwise,  are  most  important  to  the  observer  irradiance  output.  Finally,  HELEEOS 
could  be  compared  with  additional  test  data,  including  that  already  available  from 
the  December  2009  test  in  Dahlgren,  VA. 

6.1.1  Parameter  sensitivities. 

Since  the  HELEEOS  outputs  varied  so  widely  compared  with  the  measured  irra- 
diances,  it  would  be  valuable  to  conduct  a  sensitivity  study  on  the  off-axis  scattering 
model  in  order  to  determine  which  parameters  are  most  responsible  for  output  varia¬ 
tions.  Once  the  parameters  to  which  HELEEOS  is  most  sensitive  have  been  identified, 
care  should  be  taken  that  those  parameters  are  determined  as  precisely  as  possible 
while  conducting  tests  of  the  off-axis  scattering  model. 

A  high  degree  of  sensitivity  has  already  been  noted  in  the  case  of  the  off-axis  angle, 
and  this  sensitivity  is  compounded  by  its  small  magnitude  and  comparatively  large 
error  term.  This  uncertainty  may  be  reduced  somewhat  by  increasing  the  off-axis 
angle,  although  doing  so  may  make  it  more  difficult  to  observe  the  beam. 

6.1.2  Different  geometries. 

In  addition  to  increasing  the  off-axis  angle,  observations  should  be  made  looking 
further  down  the  beam  rather  than  only  a  small  distance  from  the  source.  A  small 
number  of  images  in  the  existing  data  set  show  irradiance  farther  down  the  beam, 
but  these  were  not  used  because  of  difficulty  in  accurately  determining  the  angles 
associated  with  individual  pixels  in  those  images.  With  some  modification  to  the 
calibration  software  it  should  be  possible  to  do  so. 

All  the  images  in  the  existing  data  set  showed  forward  scattering,  and  although 


71 


it  has  been  seen  that  side- scattering  may  not  be  observable,  it  should  be  possible  to 
produce  images  of  back-scattering  with  a  similar  laser  power.  Indeed,  back-scattering 
was  observed  (but  not  yet  analyzed)  in  the  recent  Dahlgren,  VA  test. 

6.1.3  Higher  power  tests. 

It  would  be  beneficial  to  test  the  off-axis  scattering  model  against  different  laser 
equipment  and  under  different  environmental  conditions.  The  July  2009  tests  were 
conducted  in  very  high  visibility  fair-weather  conditions,  with  a  narrow  wavelength 
range  and  relative  low  power  (compared  with  weaponized  lasers).  Tests  against  higher 
powered  lasers,  different  wavelengths  and  different  weather  conditions  would  be  ben¬ 
eficial.  Data  from  the  December  2009  test  in  Dahlgren,  VA,  once  they  are  available, 
will  help  fulfill  some  of  these  needs.  Although  much  of  that  test  was  also  conducted  in 
fair  weather,  the  weather  conditions  were  certainly  somewhat  different.  Furthermore, 
the  laser  in  use  was  much  more  powerful  and  operated  at  a  significantly  different 
wavelength  than  that  used  in  the  July  2009  test. 

Because  the  Dahlgren,  VA  data  were  collected  using  CCD  cameras  rather  than  an 
InGaAs  array,  and  because  of  the  different  calibration  equipment  available  on  site,  a 
different  calibration  process  must  be  used.  The  existing  calibration  software  could  be 
modified  to  accomplish  this.  The  CCD’s  used  at  Dahlgren  are  to  be  calibrated  against 
a  laser  source,  which  will  not  produce  the  constant  radiation  flux  across  the  field  of 
view  that  the  12-inch  blackbody  provided.  Thus,  an  additional  step  may  be  required 
to  measure  the  variance  in  quantum  efficiency  across  the  CCD  or,  alternatively,  it 
may  be  assumed  that  the  quantum  efficiency  is  constant  across  the  array.  In  either 
case,  the  existing  calibration  program  is  designed  specifically  to  handle  calibration 
against  a  blackbody,  and  would  likely  require  modification  to  accommodate  the  new 
calibration  process.  The  calibration  software  also  would  require  modification  in  order 
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to  be  able  to  read  the  file  format  used  to  store  the  Dahlgren  data. 

6.1.4  Improvements  to  the  calibration  software. 

The  analysis  process  revealed  some  tolerable  but  noticeable  performance  limita¬ 
tions  in  the  calibration  tool.  These  could  be  addressed  by  re-writing  the  tool  in  a 
compiled  language  such  as  C  or  C++,  or  by  implementing  caching  in  order  to  re¬ 
duce  the  number  of  floating  point  computations  required.  The  latter  would  be  more 
straightforward  since  the  graphics  features  of  the  program  use  the  matplotlib  library 
which  has  no  bindings  for  low-level  languages.  Alternatively,  the  low-level  file  i/o 
and  numerical  codes  could  be  rewritten  in  a  compiled  language  and  the  rest  left  in 
Python. 

Integrating  the  calibration  tool  with  HELEEOS  would  make  the  analysis  process 
more  convenient.  It  is  possible  to  call  Matlab  from  inside  a  Python  script,  and  this 
capability  was  used  to  call  HELEEOS  in  order  to  generate  many  of  the  numbers  used 
in  this  thesis.  However,  this  was  done  from  a  separate  script,  and  was  not  integrated 
into  the  main  calibration  tool.  Doing  so  would  be  fairly  straightforward,  depending 
on  what  subset  of  the  HELEEOS  inputs  and  outputs  were  exposed. 

Another,  potentially  more  interesting,  modification  would  be  to  add  a  motion 
tracking  capability.  This  would  provide  a  means  to  use  some  of  the  images  collected 
at  the  July  2009  test  which  could  not  be  analyzed  with  the  current  program  because 
the  camera  was  shaking  or  panning  during  the  shot.  This  would  also  improve  the 
accuracy  of  output  for  some  of  the  images  which  were  analyzed  in  spite  of  moderate 
camera  shake. 
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6.1.5  Inputting  real  weather  data  into  HELEEOS. 


The  existing  data  set  still  leaves  considerable  possibilities  for  further  study.  In 
addition  to  analyzing  some  of  the  images  omitted  from  analysis  in  this  study,  all 
the  images  could  be  re-analyzed  using  real  weather  data.  By  default,  HELEEOS 
determines  weather  conditions  from  a  climatological  database,  taking  into  account 
the  location,  season  (summer  or  winter),  time  of  day,  and  relative  humidity.  For  this 
study,  only  the  location  was  changed  from  the  default.  Future  studies  could  replace 
these  settings  with  actual  weather  conditions. 

Most  weather  data  could  simply  be  acquired  from  National  Weather  Service 
records,  but  as  discussed  previously,  the  recorded  visibilities  are  frequently  well  below 
the  actual  ones.  Instead  of  inputting  a  visibility  to  HELEEOS,  data  from  an  aerosol 
particle  counter  could  be  used  as  inputs  to  the  Mie  scattering  model.  Particle  counter 
data  are  available  for  the  current  data  set  but  were  not  used  in  the  current  study.  Af¬ 
ter  the  actual  weather  conditions  are  taken  into  account,  a  more  detailed  comparison 
between  HELEEOS  and  the  test  data  may  be  possible,  hopefully  providing  a  means 
to  identify  areas  for  improvement  in  the  HELEEOS  model. 

6.2  Closing  remarks 

The  techniques  described  here  have  proven  effective  for  measuring  off-axis  irradi- 
ances  for  comparison  with  HELEEOS.  Each  laser  scattering  image  can  be  reduced 
to  a  single  irradiance  value  which  can  be  compared  with  HELEEOS  output.  The 
present  study  found  considerable  discrepancies  between  HELEEOS  and  actual  irradi¬ 
ance  values,  but  was  unable  to  determine  whether  HELEEOS  itself  or  its  input  data 
were  most  responsible  for  the  discrepancies  in  irradiance  results.  Further  studies, 
with  these  data  or  with  others,  may  provide  a  means  to  do  so. 
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Appendix  A.  Calibration  software  user  documentation 


This  section  describes  the  calibration  software  in  more  detail  from  a  user’s  per¬ 
spective,  providing  a  general  description  of  the  interface  and  how  it  is  used. 

1.1  Dependencies 

The  software  has  been  tested  on  Ubuntu  Linux  9.10.  It  should  be  possible  to  run 
it  on  other  operating  systems  once  the  following  dependencies  are  installed: 

•  python  2.4 

•  pygtk  2.4  or  greater  with  libglade  support 

•  NumPy  1.3.0  or  greater 

•  SciPy  0.7.1  or  greater 

•  matplotlib  0.99.1  or  greater 

The  script  analyseFieldlmages.py  calls  HELEEOS  to  evaluate  predicated  irradi- 
ances;  this  script  adds  the  following  dependencies  in  order  to  interface  with  MATLAB: 

•  mlabwrap  1.1  or  greater 

•  csh  (required  by  mlabwrap) 

1.2  Using 

The  application  is  launched  by  running  calgui  .py.  This  will  display  the  applica¬ 
tion’s  main  window,  shown  in  Figure  36 
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Figure  36.  The  calibration  program’s  main  window 
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1.2.1  Calibration  scenarios. 


A  calibration  scenario  defines  all  the  blackbody  images  and  relevant  settings  used 
to  build  a  calibration,  as  well  as  field  images  to  which  the  calibration  is  to  be  applied. 
Calibration  scenarios  may  be  loaded  by  selecting  “Open”  from  the  File  menu  in  the 
main  window,  and  saved  by  selecting  “Save”  from  the  File  menu. 


1.2.2  Loading  and  editing  images. 

The  Add  Image  button  launches  a  file  chooser  dialog  in  which  one  can  pick  a 
XeniCS  XVI  file  to  load.  The  image  is  added  to  the  Blackbody  Images  table  or  the 
Field  Images  table,  depending  on  which  of  those  tabs  is  selected. 

When  an  image  is  added,  an  Image  Properties  dialog  appears,  prompting  for 
temperature,  integration  time,  and  a  frame  range  for  averaging.  The  temperature 
should  be  the  temperature  of  the  blackbody  for  blackbody  images;  it  is  not  used  for 
field  images.  All  properties  can  be  edited  later  from  the  table. 

The  Blackbody  Images  table,  shown  in  Figure  37,  shows  all  the  blackbody  images 
defined  in  the  loaded  calibration  scenario.  Table  3  lists  the  fields  in  this  table. 


Table  2.  Fields  in  the  Blackbody  Images  table 


Use 

Temperature 
BB  energy  per  pixel0 


Integration  time 

Start  frame 
Stop  frame 
Filename0 


Whether  the  image  is  used  in  computing 

the  calibration 

The  blackbody  temperature 

The  energy  in  Joules  received  by  each  pixel 

according  to  Planck  theory  for  a  blackbody 

of  the  given  temperature 

The  integration  time  set  when  taking  the 

image 

The  start  of  the  frame  range  for  averaging 
The  end  of  the  frame  range  for  averaging 
The  path  of  the  XVI  file 


“Non-editable  field 
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Figure  37.  The  table  of  blackbody  images 


The  Field  Images  table,  shown  in  Figure  38,  shows  all  the  blackbody  images 
defined  in  the  loaded  calibration  scenario.  Table  3  lists  the  fields  in  this  table. 


Table  3.  Fields  in  the  Field  Images  table 


Bg  Sub 


Integration  time 

Start  frame 
Stop  frame 
Filename® 


Whether  to  do  background  subtraction  for 
the  image  (only  works  after  the  average  has 
been  computed) 

The  integration  time  set  when  taking  the 
image 

The  start  of  the  frame  range  for  averaging 
The  end  of  the  frame  range  for  averaging 
The  path  of  the  XVI  file 


“Non-editable  field 


1.2.3  Viewing  images. 

When  an  image  is  selected  in  one  of  the  image  tables,  the  figure  pane  is  imme¬ 
diately  updated  according  the  the  view  settings.  The  data  to  be  displayed  may  be 
selected  under  the  following  options  shown  in  the  “View”  section  of  the  main  window. 
Table  4  shows  the  possible  values. 

The  frame  to  be  shown  is  selected  using  a  slider  under  the  figure  area,  or  if  the 
average  has  been  computed  the  average  may  be  displayed  instead  of  the  selected  frame 
by  selecting  the  “Averaged”  radio  button  in  the  “Display  frames”  section. 
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Figure  38.  The  table  of  field  images 


Option 

Original  Image 
Calibrated  Image 
Gain 
Offset 
Fit  variance 


Table  4.  View  types 

Data  displayed 

Raw  image  data  (possibly  averaged) 
Image  data  with  calibration  applied 
Calibration  gain  array 
Calibration  offset  array 
Errors  computed  for  each  pixel 


1.2.4  Overlays. 

An  overlay  is  an  array  of  boolean  data  which  may  be  drawn  on  top  of  the  image. 
Two  data  sets  may  be  overlaid:  The  “bad  pixel”  array  computed  from  the  calibration, 
and  the  “selected  region,”  which  is  created  according  to  the  procedure  in  Section  1.4.2. 
The  overlay  color  and  opacity  can  also  be  changed. 

When  the  “Edit  overlay”  box  is  checked,  the  overlay  may  be  edited  manually.  As 
long  as  the  figure  is  not  in  “pan/zoom”  or  “zoom  to  rectangle”  mode,  a  pixel  can  be 
added  to  the  overlay  array  with  a  left  mouse  button  click,  or  removed  with  a  right 
click.1 

1This  feature  is  currently  broken,  and  may  never  be  fixed  as  it  has  become  unnecessary  due  to 
improvements  in  the  way  the  bad  pixel  array  is  initially  computed. 
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1.3  Viewing  histograms 


By  default,  the  data  is  displayed  in  the  form  of  an  image,  but  by  selecting  “His¬ 
togram”  under  the  “Plot  type”  section,  a  histogram  may  be  displayed  instead.  When 
displaying  the  histogram,  the  number  of  bins  may  be  set  using  a  text  box.  If  the 
“Only  selection”  box  is  checked,  the  histogram  will  be  computed  using  the  pixels  of 
the  selected  region  (see  Section  1.4.2)  rather  than  the  entire  image. 

1.4  Image  settings 

The  Image  Settings  dialog  can  be  accessed  from  the  Image  Settings  option  on  the 
Settings  menu,  or  by  right-clicking  on  an  image  and  selecting  Image  Settings  from 
the  menu  that  appears.  The  options  in  the  Image  Settings  dialog  vary  depending  on 
whether  the  image  selected  is  a  blackbody  or  a  field  image. 

The  appearance  of  this  dialog  for  a  field  image  is  shown  in  Figure  39;  its  appear¬ 
ance  for  a  blackbody  image  is  shown  in  Figure  40. 

1.4.1  Image  section. 

This  section  is  available  for  both  blackbody  and  field  images.  It  displays  the 
filename  and  provides  a  field  to  set  or  change  the  image’s  integration  time. 

1.4.2  Field  Scenario  section. 

This  section  is  available  only  for  field  images.  It  provides  an  interface  for  setting 
observer  geometry  settings,  laser  wavelengths  and  power  levels,  the  frame  range,  and 
the  background  image  to  use. 
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Figure  39.  Image  Settings  dialog  for  a  field  image 
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Figure  40.  Image  Settings  dialog  for  a  blackbody  image 
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1.4.3  Optics  section. 


This  section  is  available  for  both  blackbody  and  held  images.  It  contains  settings 
related  to  the  optics  used  to  take  the  image. 
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Table  5.  Options  in  the  Image  Settings  dialog 

Image 

Integration  time  Camera  integration  time  in  //.s 

Blackbody 

Temperature  Blackbody  temperature 
Frames  to  average  Frame  range  to  average  for  use  in  calibra¬ 
tion 

Field  Scenario 

Laser  location  Location  of  the  laser  source  in  image  coor¬ 
dinates 

Distance  to  laser  Distance  from  the  laser  source  to  the  cam¬ 
era 

Angle  between  observer,  laser,  and  target  Angle  between  laser  beam  and  observer 

path  (angle  B  in  Figure  7) 

FOV  starting  distance  Start  distance  in  meters  along  laser  beam 

path  (will  be  converted  to  image  space) 
FOV  stopping  distance  Stop  distance  in  meters  along  laser  beam 

path  (will  be  converted  to  image  space) 
Apparent  beam  direction  Direction  beam  points  in  image  space 

(clockwise  where  0°  points  to  the  right) 
Selection  extent  How  many  degrees  off  of  the  beam  direc¬ 
tion  a  pixel  must  be  before  it  is  considered 
outside  of  the  selection 

Laser  beams  Power  and  wavelength  of  each  beam  visible 
in  the  selection  region 

Scenario  frames  Range  of  frames  to  use  when  computing 
the  irradiance  of  the  image 

Background  frames  Range  of  frames  to  average  for  use  as  a 

background 

Background  image  Filename  of  the  hie  to  use  for  the  back¬ 
ground.  Defaults  to  the  current  image. 

Optics 

Use  optics  settings  from  calibration  Whether  the  optics  settings  (aperture,  ob¬ 
scuration,  FOV)  are  the  same  as  were  used 
in  the  calibration  (other  fields  in  this  sec¬ 
tion  are  disabled  if  this  is  checked) 
Aperture  Aperture  diameter  of  the  receiving  optic 
FOV  per  pixel  Field  of  view  in  radians  of  a  single  pixel 
Fractional  obscuration  Fraction  of  the  aperture  that  is  obscured 

by  a  secondary  mirror  or  other  objects 


84 


Appendix  B.  HELEEOS  inputs 


The  HELEEOS  settings  given  here  were  used  by  default  for  all  HELEEOS  runs 
in  the  project.  Observer  and  laser  settings  were  overridden  from  these  defaults  based 
on  the  conditions  under  which  each  field  image  was  taken.  Specifically,  the  obserever 
“Relative  Azimuth  from  Platform”  setting  was  replaced  by  (90  —  0obs)  where  0oiJS  is  the 
observer’s  off-axis  angle,  and  the  observer  “Initial  Distance  from  Platform”  setting 
was  replaced  by  the  actual  distance  from  the  source  laser  to  the  observer. 

Table  6.  HELEEOS  inputs 


Scenario:  Location  tab 


Latitude 

39.83 

Longitude 

-84.05 

ExPERT  Site 

DAYTON /WRIGHT-PATTE 

Season 

Summer 

Time  of  Day 

00:00-03:00 

ExPERT  percentile 

50%  -  Average 

Enable  Advanced  Atmosphere  Settings 

Checked 

Use  Fast  Atmospheric  Calculations 

Unchecked 

Select  Atmosphere 

ExPERT  Atmospheres 

Enable  User  Specified  Ground  Visibility 

Checked 

Aerosol  Visibility  (km) 

100 

Select  Molecular  Effects  Model 

Calculations 

Select  Aerosol  Effects  Model 

GADS  Worldwide  Aerosols 

Enable  User  Specified  Boundary  Layer  Height 

Unchecked 
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Scenario:  Atmosphere  tab 

Turbulence 

Turbulence  Multiplier 

Select  Wind  Model 

Enable  Clouds  and  Rain 

HV  5/7 

1 

Climatological  Wind  Profile 

Unchecked 

Platform:  Geometry  tab 

Initial  Altitude 

2 

Initial  Relative  Azimuth  to  next  Object 

90 

Initial  Distance  to  next  Object 

2150 

Horizontal  Velocity  (m/s) 

0 

Velocity  Heading  (0  is  North) 

320 

Vertical  Velocity  (m/s) 

0 

Vertical  Acceleration  (m/s  "2) 

0 

Platform:  Laser  tab 

Laser  Propagation  Model 

Top  Hat 

Laser  Type 

Continuous  Wave  (CW) 

Laser  Wavelength 

User  Wavelength 

User  Wavelength 

1.533e-6 

Power  (W) 

0.6 
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Platform:  Optics  tab 

Top  Hat  Focus  Settings 

Spot  Size 

Beam  Width 

l/e"2  Definition  of  Beam  Width 

Spot  Diameter  (m) 

1.6 

Relative  Obscuration 

0 

Exit  Aperture  Diameter  (m) 

0.004 

Wavefront  Error  (waves) 

0 

Beam  Quality 

1 

Total  System  RMS  Jitter  (rad) 

0 

Optics  AO 

Default 

Aero-Optic  Effect 

No  Aero  Optic  Model 

Platform:  Lethality  tab 

Shape  of  the  susceptible  area  of  the  target 

Rectangular 

Susceptible  Target  Width  (m) 

0.05 

Distribution  Type 

Normal  Distribution 

Target  Damage  Threshold  (J/nC2) 

5e7 

Target  Damage  Threshold  (J/nC2) 

le7 

Observer:  Geometry  tab 

Enable  Observer 

Checked 

Initial  Altitude  (m) 

2 

Initial  Distance  from  Platform  (m) 

2160 

Relative  Azimuth  from  Platform 

88.78 

Horizontal  Velocity  (m/s) 

0 

Vertical  Velocity  (m/s) 

0 

Vertical  Acceleration  (m/s  "2) 

0 
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Observer:  Optics  tab 

Beam  to  be  Observed 

Platform  to  Target 

Focusing 

Linear  Focusing 

FOV  Starting  Dist  (m) 

300 

FOV  Ending  Dist  (m) 

750 

Jitter  (rad) 

0 

Relative  Obscuration 

0 

Input  Aperture  Diameter  (m) 

0.101 

Observer  AO 

Default 

Relay  1:  Geometry  tab 

Enable  Relay  1 

Unchecked 

Relay  2:  Geometry  tab 

Enable  Relay  2 

Unchecked 

Edit  Libraries:  Atmosphere  tab 

Internet  Database 

Please  Select  ... 

Edit  Libraries:  Advanced  AO  tab 

Type  of  AO  System 

HWFS_DM 

Error  Rejection  Bandwidth  (Hz) 

0 

Latency  of  Transmit  Loop  (s) 

0 

Deformable  Mirror  Actuator  Spacing  (m) 

0 

Diameter  of  Beacon  from  Transmitter  (m) 

0 

Type  of  Track  System 

FPA_FSM 

Error  Rejection  Bandwidth  (Hz) 

0 

Latency  of  Transmit  Loop  (s) 

0 

Appendix  C.  Table  of  field  images  with  corresponding  data 


For  each  image  used,  the  relevant  configuration  information  is  given,  along  with 
the  irradiances  computed  by  the  Digital  Camera  Calibration  Tool  (DCCT)  and  by 
HELEEOS.  Some  image  files  are  listed  more  than  once  because  several  different  frame 
ranges  were  used,  or  because  multiple  selection  areas  were  used  due  to  lasers  pointing 
in  different  directions. 

Two  optical  setups  were  used,  referred  to  by  the  names  “Small  telescope”  and 
“Large  telescope.”  Their  properties  are  listed  in  Table  7. 

Table  7.  Properties  of  the  two  optical  setups. 

Name  Aperture  (cm)  FOV  per  pixel(rad)  Relative  obscuration 
Small  telescope  10.10  4.10628019324  x  10"5  OR 

Large  telescope  31.75  8.81752780784  x  10"6  0.130 

The  distance  from  the  laser  to  the  observer  is  given;  this  value  was  assigned  in 
the  HELEEOS  observer  settings  accordingly.  The  off-axis  angle  was  subtracted  from 
90°  and  passed  to  HELEEOS  as  the  “Relative  Azimuth  from  Platform.”  For  instance, 
where  the  off-axis  angle  was  0.9°,  the  “Relative  Azimuth  from  Platform”  parameter 
was  set  to  90  —  0.9  =  89.1°.  The  FOV  range  values  given  were  applied  to  the  “FOV 
Starting  Dist”  and  “FOV  Ending  Dist”  parameters  in  HELEEOS.  One  HELEEOS 
run  was  performed  for  each  laser  in  each  image,  and  the  results  were  added  together 
to  obtain  the  HELEEOS  irradiances  given  here. 

The  uncertainties  for  the  DCCT  output  were  computed  as  the  standard  deviation 
of  the  irradiances  for  all  the  frames. 
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Table  8.  Field  irradiance  and  corresponding  HELEEOS  output  (red  areas  in  the  images 
denote  the  region  over  which  summation  was  performed) 


2009-07-15/ 0004. xvi  (frames  1-2876) 

Optics:  Large  telescope 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

1.5887 /im 

5.50  W 

Lasers: 

1.5953/7777 

5.70  W 

1.5326/7777 

5.50  W 

DCCT  irradiance: 

1.617  x  lO"05 

±  3.07  x  10-°7  W/m2 

HELEEOS  irradiance: 

8.874  x  IQ”07 

W/m 2 

2009-07-15/ 0004.xvi  (frames  3245-6455) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

Lasers: 

1.5887/7777 

5.50  W 

1.5953/7777 

5.70  W 

DCCT  irradiance: 

1.404  x  10-°5 

±  4.95  x  10-°7  W/m2 

HELEEOS  irradiance: 

8.874  x  10“07 

W/m 2 
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2009-07-15/ 0004.xvi 

(frames  6490-9425) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Lasers: 

Wavelength  Power 

1.5887 /im  5.50  W 

DCCT  irradiance: 

4.156  x  10-°6  ±  1.77  x  10"07 

W/m 2 

HELEEOS  irradiance: 

4.386  x  IQ”07  W/m2 

2009-07-15/ 0005.xvi  (frames  4118-13056) 

Optics: 

Small  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

299  -  1200  m 

Lasers: 

Wavelength  Power 

1.5887 5.50  W 

DCCT  irradiance: 

9.421  x  10“08  ±  9.28  x  10“09 

W/m 2 

HELEEOS  irradiance: 

2.282  x  10"06  W/m2 
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le-11 


2009-07-15/ 0006.xvi 

(frames  11822-22138) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength  Power 

Lasers: 

1.5887 iim  0.60  W 

1.5953 pm  0.60  W 

DCCT  irradiance: 

9.113  x  10“08  ±  4.57  x  10“09  W/m2 

HELEEOS  irradiance: 

9.509  x  10"08  W/m2 

2009-07-15/ 0008.xvi 

(frames  7986-16537) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength  Power 

Lasers: 

1.5857^™  0.60  W 

DCCT  irradiance: 

4.301  x  lO"06  ±  1.66  x  10-°7  W/m2 

HELEEOS  irradiance: 

4.795  x  10“08  W/m2 
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le-11 


le-ll 


2009-07-15/1109.xvi  (frames  1-9273) 

Optics:  Large  telescope 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

Lasers: 

1.5326 /im 

0.60  W 

1.5857 /im 

0.60  W 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

3.464  x  10-°8 

±  7.38  x  10"09  W/m2 

HELEEOS  irradiance: 

1.444  x  10“07 

W/m2 

2009-07-15/2243.xvi  (frames  1-5588) 

Optics:  Large  telescope 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  299  m 

Wavelength 

Power 

Lasers: 

1.5326 /im 

0.60  W 

1.5857 pm 

0.60  W 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

4.345  x  10"08 

±2.13  x  lO"09  W/m2 

HELEEOS  irradiance: 

9.248  x  lO”08 

W/m2 
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le-11 


2009-07-15/2250.xvi  (frames  1-6962) 

Optics:  Large  telescope 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

Lasers: 

1.5326 /im 

0.60  W 

1.5857 /im 

0.60  W 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

6.822  x  10"08 

±  2.54  x  10-°9  W/m2 

HELEEOS  irradiance: 

1.444  x  10“07 

W/m2 

2009-07-15/2300.xvi  (frames  1-2408) 

Optics:  Large  telescope 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

Lasers: 

1.5857 pm 

0.60  W 

1.5953 pm 

0.60  W 

1.5326 /im 

0.60  W 

DCCT  irradiance: 

6.935  x  10"06 

±3.16  x  lO"07  W/m2 

HELEEOS  irradiance: 

9.519  x  lO”08 

W/m2 
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3.0 


1.8 


1.2 


2009-07-15/2300.xvi 

(frames  2583-2791) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength  Power 

Lasers: 

1.5857 iim  0.60  W 

1.5953 pm  0.60  W 

DCCT  irradiance: 

6.682  x  10“06  ±  2.72  x  10“07  W/m 2 

HELEEOS  irradiance: 

9.519  x  10-°8  W/m2 

2009-07-15/2300.xvi 

(frames  2837-3082) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength  Power 

Lasers: 

1.5857^™  0.60  W 

DCCT  irradiance: 

2.096  x  10-°6  ±  7.37  x  lO"08  W/m2 

HELEEOS  irradiance: 

4.795  x  10“08  W/m2 
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2009-07-15/2306.xvi  (frames  4571-21983) 


Optics:  Large  telescope 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

Lasers: 

1.5326 /im 

0.60  W 

1.5857 /im 

0.60  W 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

9.283  x  10"08 

±  4.70  x  10-°9  W/m2 

HELEEOS  irradiance: 

1.444  x  10“07 

W/m 2 

2009-07-15/2310.xvi  (frames  1-2050) 

Optics:  Large  telescope 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

Lasers: 

1.5326 /im 

0.60  W 

1.5857 pm 

0.60  W 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

1.230  x  10-°7 

±  4.98  x  lO"09  W/m2 

HELEEOS  irradiance: 

1.444  x  IQ”07 

W/m2 
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2009-07-15/2315-SmallTele.xvi  (frames  1-2169) 

Optics:  Small  telescope 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

299  -  1199  m 

Wavelength 

Power 

Lasers: 

1.5326 /im 

0.60  W 

1.5857 /im 

0.60  W 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

1.081  x  10-°7 

±  4.46  x  10-°9  W/m2 

HELEEOS  irradiance: 

7.515  x  10“07 

W/m 2 

2009-07-15/2318.xvi  (frames  1-4107) 

Optics:  Large  telescope 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

Lasers: 

1.5326 /im 

0.60  W 

1.5857 pm 

0.60  W 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

5.174  x  10"08 

±  2.56  x  lO"09  W/m2 

HELEEOS  irradiance: 

1.444  x  IQ”07 

W/m2 
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2009-07-15/2331.xvi  (frames  1-1663) 

Optics:  Large  telescope 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

Lasers: 

1.5326 /im 

0.60  W 

1.5857 /im 

0.60  W 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

7.887  x  10"08 

±  4.49  x  10-°9  W/m 2 

HELEEOS  irradiance: 

1.444  x  10“07 

W/m2 

2009-07-15/2334. xvi  (frames  1-1434) 

Optics:  Large  telescope 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

Lasers: 

1.5857 pm 

0.60  W 

1.5326/um 

0.60  W 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

7.433  x  10"08 

±  1.64  x  lO"09  W/m2 

HELEEOS  irradiance: 

4.795  x  lO”08 

W/m2 
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2009-07-15/2334.xvi 

(frames  1510-1612) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength  Power 

Lasers: 

1.5857 iim  0.60  W 

1.5953 pm  0.60  W 

DCCT  irradiance: 

6.837  x  lO”08  ±  1.47  x  10“09  W/m2 

HELEEOS  irradiance: 

4.795  x  10"08  W/m2 

2009-07-15/2334.xvi 

(frames  1623-1910) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength  Power 

Lasers: 

1.5857^™  0.60  W 

DCCT  irradiance: 

2.199  x  lO"08  ±  6.52  x  lO"10  W/m2 

HELEEOS  irradiance: 

4.795  x  10“08  W/m2 
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2009-07-15/2337.xvi  (frames  3966-6004) 


Optics: 


Large  telescope 


50  100  150  200  250  300 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

Lasers: 

1.5857 /im 

0.60  W 

DCCT  irradiance: 

6.077  x  10"09 

±  3.73  x  10"10  W/m2 

HELEEOS  irradiance: 

4.795  x  lO”08 

W/m2 

2009-07-15/2340.xvi  (frames  1-3224) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

1.5326/um 

0.60  W 

Lasers: 

1.5857 fj,m 

0.60  W 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

7.711  x  10“08 

±  2.20  x  10“09  W/m2 

HELEEOS  irradiance: 

1.444  x  lO"07 

W/m 2 

50  100  150  200  250  300 
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2009-07-15/2351.xvi  (frames  1-2862) 

Optics:  Large  telescope 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

Lasers: 

1.5326 /im 

0.60  W 

1.5857 /im 

0.60  W 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

8.053  x  10-°8 

±  4.07  x  10"09  W/m2 

HELEEOS  irradiance: 

1.444  x  10“07 

W/m2 

2009-07-15/2354.xvi  (frames  1-3690) 

Optics:  Large  telescope 


Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

Lasers: 

1.5857 pm 

0.60  W 

1.5326/um 

0.60  W 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

8.093  x  10"08 

±  2.61  x  lO"09  W/m2 

HELEEOS  irradiance: 

1.444  x  IQ”07 

W/m2 
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2009-07-15/2354.xvi 

(frames  3717-3874) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength  Power 

Lasers: 

1.5857 iim  0.60  W 

1.5953 pm  0.60  W 

DCCT  irradiance: 

7.570  x  10“08  ±  1.98  x  10“09  W/m2 

HELEEOS  irradiance: 

1.444  x  10-°7  W/m 2 

2009-07-15/2354.xvi 

(frames  3896-4014) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength  Power 

Lasers: 

1.5857^™  0.60  W 

DCCT  irradiance: 

2.534  x  lO"08  ±  7.91  x  lO"10  W/m 2 

HELEEOS  irradiance: 

1.444  x  10“07  W/m2 
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2009-07-15/2356.xvi 

(frames  5536-12133) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

Lasers: 

1.5953 pm 

5.40  W 

DCCT  irradiance: 

5.133  x  10-°6 

±  2.86  x  lO"07  W/m2 

HELEEOS  irradiance: 

4.252  x  IQ”07 

W/m2 

2009-07-15/2356.xvi  (frames  12452-13903) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

1.5326/um 

5.70  W 

Lasers: 

1.5857 fj,m 

5.60  W 

1.5953 pm 

5.40  W 

DCCT  irradiance: 

3.904  x  lO”05 

±  1.31  x  10“06  W/m2 

HELEEOS  irradiance: 

1.340  x  10"06 

W/m 2 
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2009-07-15/2358.xvi 

(frames  1-21247) 

Optics: 

Large  telescope 

Distance: 

3015  m 

Off-axis  angle: 

0.900° 

FOV  range: 

100  -  400  m 

Wavelength 

Power 

1.5326 pm 

5.50  W 

Lasers: 

1.5857 /im 

5.50  W 

1.5953 pm 

5.70  W 

DCCT  irradiance: 

5.828  x  10-°7 

±  6.60  x  10-°8  W/m2 

HELEEOS  irradiance: 

1.340  x  10“06 

W/m2 

2009-07-17 / 0700.xvi  (frames  4185-7401) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  749  m 

Wavelength 

Power 

Lasers: 

1.5326/rm 

0.60  W 

DCCT  irradiance: 

3.601  x  10“08 

±  5.65  x  10“09  W/m2 

HELEEOS  irradiance: 

2.627  x  10"07 

W/m2 
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2009-07-17 / 0708.xvi 

(frames  1-5178) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  749  m 

Wavelength  Power 

Lasers: 

1.5326 pm  0.60  W 

DCCT  irradiance: 

2.974  x  10-°8  ±  2.49  x  10"09  W/m 2 

HELEEOS  irradiance: 

2.626  x  IQ”07  W/m2 

2009-07-17 / 0712. xvi  (frames  1-6274) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  749  m 

Wavelength  Power 

Lasers: 

1.5887 0.60  W 

DCCT  irradiance: 

2.494  x  10“08  ±  1.86  x  10“°9  W/m2 

HELEEOS  irradiance: 

2.521  x  IQ"07  W/m2 
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2009-07-17 / 0725.xvi 

(frames  3685-11954) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  749  m 

Lasers: 

Wavelength  Power 

1.5326 pm  0.60  W 

DCCT  irradiance: 

2.184  x  10-°8  ±  2.52  x  10"09 

W/m2 

HELEEOS  irradiance: 

2.626  x  IQ”07  W/m2 

2009-07-17 / 0728.xvi  (frames  1-4755) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  750  m 

Lasers: 

Wavelength  Power 

1.5887 0.60  W 

DCCT  irradiance: 

1.607  x  10“08  ±  1.66  x  10“09 

W/m 2 

HELEEOS  irradiance: 

2.522  x  IQ"07  W/m,2 
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2009-07-17 / 0730.xvi 

(frames  1-8817) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  749  m 

Lasers: 

Wavelength  Power 

1.5326 pm  0.60  W 

DCCT  irradiance: 

2.075  x  10-°8  ±  1.93  x  10"09 

W/m2 

HELEEOS  irradiance: 

2.626  x  IQ”07  W/m2 

2009-07-17 / 0730.xvi  (frames  2049-5952) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  749  m 

Lasers: 

Wavelength  Power 

1.5887 0.60  W 

DCCT  irradiance: 

2.532  x  10“08  ±  9.00  x  10“°9 

W/m2 

HELEEOS  irradiance: 

2.520  x  IQ"07  W/m,2 
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2009-07-17 / 0730.xvi 

(frames  5981-8817) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  749  m 

Wavelength  Power 

Lasers: 

1.5887 iim  0.60  W 

1.5953 pm  0.60  W 

DCCT  irradiance: 

2.532  x  10“08  ±  9.01  x  10“°9  W/m 2 

HELEEOS  irradiance: 

2.520  x  10-°7  W/m2 

2009-07-17 / 0734.xvi 

(frames  1-3147) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  750  m 

Wavelength  Power 

Lasers: 

1.5326 nm  0.60  W 

DCCT  irradiance: 

2.300  x  10"08  ±  8.71  x  10"09  W/m2 

HELEEOS  irradiance: 

2.628  x  10“07  W/m2 
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2009-07-17 / 0734.xvi 

(frames  1-3147) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  750  m 

Wavelength  Power 

Lasers: 

1.5887 /im  0.60  W 

1.5953 pm  0.60  W 

DCCT  irradiance: 

3.844  x  10“08  ±  8.64  x  10“°9  W/m2 

HELEEOS  irradiance: 

2.522  x  10-°7  W/m2 

2009-07-17 / 0736.xvi 

(frames  1-5702) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  749  m 

Lasers: 

Wavelength  Power 

1.5326 urn  0.60  W 

DCCT  irradiance:  1.967  x  10  08  ±  1.81  x  10  09  W/m 2 

HELEEOS  irradiance:  2.627  x  10-07  W/m 2 
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2009-07-17 / 0736.xvi 

(frames  1853-5702) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  749  m 

Lasers: 

Wavelength  Power 

1.5953 pm  0.60  W 

DCCT  irradiance: 

1.624  x  10-°8  ±  1.22  x  10-°9 

W/m2 

HELEEOS  irradiance: 

2.489  x  IQ”07  W/m2 

2009-07-17 / 0736.xvi  (frames  1936-5702) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  749  m 

Lasers: 

Wavelength  Power 

1.5887 0.60  W 

DCCT  irradiance: 

8.832  x  10“09  ±  1.73  x  10“°9 

W/m2 

HELEEOS  irradiance: 

2.521  x  IQ"07  W/m,2 

110 


le-10 


le-10 


0  50  100  150  200  250  300 


= 

1.6 


0.8 


2009-07-17 / 0737.xvi 

(frames  1-3224) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  750  m 

Wavelength  Power 

Lasers: 

1.5326 pm  0.60  W 

DCCT  irradiance: 

2.101  x  10-°8  ±  1.67  x  10"09  W/m2 

HELEEOS  irradiance: 

2.627  x  10“07  W/m2 

2009-07-17 / 0737.xvi  (frames  1-3224) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  750  m 

Wavelength  Power 

Lasers: 

1.5953/um  0.60  W 

DCCT  irradiance: 

4.868  x  10“09  ±  1.67  x  10“09  W/m2 

HELEEOS  irradiance: 

2.490  x  IQ"07  W/m2 
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2009-07-17 / 0737.xvi 

(frames  1-3224) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  750  m 

Lasers: 

Wavelength  Power 

1.5887 /im  0.60  W 

DCCT  irradiance: 

1.145  x  10-°8  ±  1.31  x  10-°9 

W/m2 

HELEEOS  irradiance: 

2.521  x  IQ”07  W/m2 

2009-07-17 / 0738.xvi  (frames  1-6334) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

300  -  749  m 

Lasers: 

Wavelength  Power 

1.5887 0.60  W 

DCCT  irradiance: 

8.973  x  10“09  ±  2.59  x  10“°9 

W/m2 

HELEEOS  irradiance: 

1.831  x  IQ"07  W/m,2 
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2009-07-17 / 0738.xvi 

(frames  1-9193) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  749  m 

Lasers: 

Wavelength  Power 

1.5953 pm  0.60  W 

DCCT  irradiance: 

7.416  x  10"09  ±  1.59  x  10"09 

W/m2 

HELEEOS  irradiance: 

2.453  x  IQ”07  W/m2 

2009-07-17 / 0738.xvi  (frames  1-6060) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  749  m 

Lasers: 

Wavelength  Power 

1.5326 pm  0.60  W 

DCCT  irradiance: 

1.774  x  10“08  ±  1.45  x  10“°9 

W/m2 

HELEEOS  irradiance: 

2.627  x  10"07  W/m2 
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2009-07-17 / 0742.xvi  (frames  1-889) 

Optics:  Small  telescope 


Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  750  m 

Wavelength 

Power 

Lasers: 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

1.206  x  10-°6 

±  9.62  x  10-°8  W/m2 

HELEEOS  irradiance: 

5.150  x  IQ”07 

W/m2 

2009-07-17 / 0742.xvi  (frames  1-889) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  750  m 

Wavelength 

Power 

Lasers: 

1.5326/um 

0.60  W 

1.5887 fxm 

0.60  W 

DCCT  irradiance: 

1.229  x  10-°6 

±  2.83  x  10-°8  W/m2 

HELEEOS  irradiance: 

5.150  x  IQ”07 

W/m 2 
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2009-07-17 / 0806.xvi 

Optics: 

Distance: 

Off-axis  angle: 

FOV  range: 


(frames  1-3224) 

Small  telescope 
2160  m 
1.220° 

99  -  750  m 


0  50  100  150  200  250  300  ■ -2.0 


Lasers: 

Wavelength  Power 

1.5953 pm  0.60  W 

DCCT  irradiance: 

-2.386  x  10"08  ±  2.85  x  10"08  W/m2 

HELEEOS  irradiance: 

2.491  x  IQ”07  W/m2 

2009-07-17 / 0806.xvi  (frames  1-3224) 

Optics: 

Small  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

99  -  750  m 

Wavelength  Power 

Lasers: 

1.5326/um  0.60  W 

1.5887 urn  0.60  W 

DCCT  irradiance: 

6.042  x  10-°8  ±  2.70  x  10"08  W/m 2 

HELEEOS  irradiance: 

5.150  x  IQ”07  W/m 2 

50  100  150  200  250  300 
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2009-07-17/0811.xvi  (frames  1-3224) 

Optics:  Large  telescope 


Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

24  -  199  m 

Wavelength 

Power 

Lasers: 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

1.106  x  10-°8 

±  8.02  x  10"10  W/m2 

HELEEOS  irradiance: 

4.567  x  10“08 

W/m 2 

2009-07-17/0811.xvi  (frames  1-3224) 

Optics: 

Large  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

50  -  225  m 

Wavelength 

Power 

Lasers: 

1.5887 fxm 

0.60  W 

1.5326 pm 

0.60  W 

DCCT  irradiance: 

1.554  x  10-°8 

±  9.20  x  10"10  W/m2 

HELEEOS  irradiance: 

9.773  x  lO”08 

W/m 2 
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2009-07-17 / 0815.xvi 

(frames  1-6585) 

Optics: 

Large  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

24  -  219  m 

Wavelength  Power 

Lasers: 

1.5887 iim  0.60  W 

1.5326 pm  0.60  W 

DCCT  irradiance: 

3.578  x  10“08  ±  7.17  x  10“°9  W/m2 

HELEEOS  irradiance: 

1.082  x  10-°7  W/m2 

2009-07-17 / 0815.xvi 

(frames  1-6585) 

Optics: 

Large  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

24  -  220  m 

Wavelength  Power 

Lasers: 

1.5953^7/1  0.60  W 

DCCT  irradiance: 

2.569  x  lO"08  ±  5.29  x  lO"09  W/m2 

HELEEOS  irradiance: 

5.239  x  10“08  W/m2 
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2009-07-17 / 0817.xvi 

(frames  1-20234) 

Optics: 

Large  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

49  -  225  m 

Wavelength  Power 

Lasers: 

1.5887 iim  0.60  W 

1.5326 pm  0.60  W 

DCCT  irradiance: 

1.984  x  10“08  ±  8.37  x  10“09  W/m2 

HELEEOS  irradiance: 

9.843  x  10"08  W/m2 

2009-07-17 / 0817.xvi 

(frames  1-20234) 

Optics: 

Large  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

24  -  220  m 

Wavelength  Power 

Lasers: 

1.5953 nm  0.60  W 

DCCT  irradiance: 

1.955  x  10"08  ±  6.82  x  10"09  W/m2 

HELEEOS  irradiance: 

5.239  x  10“08  W/m2 
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le-10 


le-10 


1.2 

1.0 

0.8 


0.0 


2009-07-17 / 0819.xvi 

(frames  1-3224) 

Optics: 

Large  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

49  -  224  m 

Wavelength  Power 

Lasers: 

1.5326 pm  0.60  W 

1.5887 /im  0.60  W 

DCCT  irradiance: 

7.497  x  10“08  ±  3.29  x  10“°9  W/m2 

HELEEOS  irradiance: 

9.829  x  10"08  W/m2 

2009-07-17 / 0819.xvi 

(frames  1-3224) 

Optics: 

Large  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

24  -  220  m 

Wavelength  Power 

Lasers: 

1.5953^7/1  0.60  W 

DCCT  irradiance: 

6.986  x  lO"08  ±  2.86  x  lO"09  W/m2 

HELEEOS  irradiance: 

5.147  x  10“08  W/m2 

119 


2009-07-17 / 0820.xvi  (frames  1-3224) 

Optics:  Large  telescope 


Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

24  -  220  m 

Wavelength 

Power 

Lasers: 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

2.657  x  10"08 

±  9.56  x  lO"10  W/m2 

HELEEOS  irradiance: 

5.227  x  10“08 

W/m 2 

2009-07-17 / 0820.xvi  (frames  1-3224) 

Optics: 

Large  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

49  -  224  m 

Wavelength 

Power 

Lasers: 

1.5326/um 

0.60  W 

1.5887 fxm 

0.60  W 

DCCT  irradiance: 

2.718  x  10-°8 

±  9.34  x  lO"10  W/m2 

HELEEOS  irradiance: 

9.829  x  lO”08 

W/m 2 
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2009-07-17 / 0843.xvi  (frames  1-2232) 

Optics:  Large  telescope 


le-ll 


3.2 


Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

24  -  199  m 

Wavelength 

Power 

Lasers: 

1.5953 pm 

0.60  W 

DCCT  irradiance: 

3.917  x  10"09 

±  1.01  x  lO"09  W/m2 

HELEEOS  irradiance: 

4.646  x  lO”08 

W/m2 

2009-07-17 / 0843.xvi  (frames  1-1594) 

Optics: 

Large  telescope 

Distance: 

2160  m 

Off-axis  angle: 

1.220° 

FOV  range: 

49  -  224  m 

Wavelength 

Power 

Lasers: 

1.5887 [im 

0.60  W 

1.5326/um 

0.60  W 

DCCT  irradiance: 

5.660  x  10"09 

±  4.25  x  10"10  W/m2 

HELEEOS  irradiance: 

9.829  x  lO”08 

W/m 2 
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